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Abstract 

The Schwarzschild and Reissner-Nordstr0m solutions to Einstein's equations describe space- times 
which contain spherically symmetric black holes. We consider solutions to the linear wave equation in 
the exterior of a fixed black hole space- time of this type. We show that for solutions with initial data 
which decay at infinity and at the bifurcation sphere, a weighted norm in space decays like . This 
weight vanishes at the event horizon, but not at infinite. 
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1 Introduction 



Black holes are very important objects in Relativity, but very little is known about their dynamics and 
interaction. Even the question of stability under small perturbations is a challenging, open problem. Any 
solution to such questions will require an understanding of the interaction between black holes and grav- 
itational radiation. The structure of the vacuum Einstein equations, which govern these dynamics in the 
absence of other matter, make it impossible to consider only radial perturbations, so that more general 
perturbations must be considered immediately. We hope that the study of linear, uncoupled waves outside 
the black hole will help provide an understanding of these problems. 

Even the stability of the empty space was a very difficult problem. It was solved originally by Christodoulou 
and Klainerman 7 , and Lindblad and Rodnianski have since developed a simpler proof 23 . Decay estimates 
for solutions of non linear wave equations played an important role in both proofs. 

In Relativity, the structure of space- time is determined by a Lorentzian pseudo- metric, g, which, in the 
absence of matter, satisfies the Einstein equations, 

R^.-^Rg^.^O. (1.1) 

The Schwarzschild and Reissner-Nordstr0m solutions are singular solutions to the Einstein equations 
which describe space- times containing a spherically symmetric black hole. The Schwarzschild solution is 
a special case of the Reissner-Nordstr0m solutions. We will restrict our attention to the exterior region, 
outside the black hole. We discuss the geometry of the exterior region further in subsection l2.1l Because the 
structure of the Einstein equations makes it impossible for these to have spherically symmetric perturbations, 
it is expected that the study of black- hole stability will require analysis of the more general, non spherically 
symmetric, rotating, Kerr-Newman black holes. The linearisation of the Einstein equations around the Kerr 
solution has been shown to have no unstable modes 

In 1957, Regge and Wheeler investigated the linear stability of the Schwarzschild black- hole, and were able 
to reduce the problem to the study of a second order, scalar equation 27 . After appropriate transformations, 
this equation and the geometrically defined, scalar wave equation differ only by a multiple of the potential 
term appearing in each. Using these transformations, the exterior region of the black- hole can be decomposed 
into the product of time and a three dimensional space. 

Because of the intrinsic interest in the geometric wave equation and because of its possible applications 
to the study of black hole stability, we consider the wave equation 

□m=0, (1.2) 

u(l) =uo, (1.3) 

^(1) -ui. (1.4) 

For this equation, we prove the following result. 

Theorem 18.211 // u is a solution to the wave equation H2.21|l , and u = ru, 

WPiuhom <t'^Ci\\uo\\l + E[uo,ui]+Ec[uo,ui] + \\L'u\\^ + E[L'uo, L'ui])^ , (1.5) 
\\{pl + l)^u\\LHm) <t'~'C{\\uo\\l + E[uo,ui] + Ec[uo, m] + \\L'u\\^ + E[L'uo, L'ui])^ (1.6) 

where and Ec[v,w] are the energy and conformal charge which are defined in section\^ is the 

Regge-Wheeler radial co-ordinate defined in section\^ and L is an angular derivative operator defined in 
section^ The measures are defined in section\^ 

Previously, we were able to apply our techniques to prove results for both the wave equation on the 
Schwarzschild manifold and the Regge-Wheeler equation, and expect the same to be true here^E]. 

This decay rate is slower than the rate found in Euclidean space, which is |^. This is consistent with 
current results for other equations on the Schwarzschild solution, which only prove slower rates of decay than 
those found in R^+^. The L°° norm has been shown to decay like t~e- for both the massive Dirac equation 
[T7j and the massive Klein- Gordon equation (2l4 . For the Schrodinger equation, the L°° has been shown to 
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decay like All these results hold only for radial initial data, or data which is the sum of finitely 

many spherical harmonics. 

The rigorous study of scattering for linear waves on black hole backgrounds was begun by Dimock and 
Kay who proved existence and completeness of wave operators for the linear wave equation ^21 and the Klein- 
Gordon equationjJJ. Asymptotic completeness and global existence for the cubic semi linear wave equation 
have been proven for a wide class of black hole manifolds, including the Kerr- Newman solutions^l, ,26,. 
DeBievre, Hislop, and Sigal have proven asymptotic completeness on a large class of non compact manifolds, 
including the Reissner-Nordstr0m manifolds ITT]. However, none of these methods give decay estimates even 
in the linear case. For sufficiently super critical Reissner-Nordstr0m black holes, Strichartz estimates have 
also been proven[B]. Little is known about fields coupled to the Einstein equations, but for a spherically 
symmetric, scalar field, Dafermos and Rodnianski have shown the field decays along the event horizon 
according to an inverse cubic. Price law They also proved a weaker rate of decay for decoupled, radial, 
semi linear wave equations in the exterior of Reissner-Nordstr0m solutions jJUj . 

We start by introducing an analogue of the conformal charge used by Ginibre and Velo ^Hl- However, 
we are prevented from completing the argument which holds for the wave equation in Euclidean space by 
the presence of a closed geodesic surface in the three dimensional space which describes the exterior region 
of the black hole. The absence of such geodesies is a non- trapping condition which is commonly imposed in 
scattering theory. The presence of this geodesic surface and the gravitational lensing they cause is already 
known [51 15^ . This surface is called the photon sphere. 

To over come this obstacle, it is sufficient to prove estimates on the angular derivative of the solution in 
the region near the closed geodesies. We prove: 

Theorem 18.21)1 If e > 0, then if u is a solution to the wave equation (|2.21|) and u = ru 

\\L'-ixc.u\\^dt < C{\\uo\\l2+E[uo,ui]), (1.7) 



where L acts like one angular derivative and Xa is a function with compact support near the closed geodesies. 

This is sufficient to prove theorem 18.211 The loss of L'^ is responsible for the additional factors of 
appearing in theorem 18.211 

Estimates on the angular derivative have already been used to prove Strichartz and point wise in time 
LP estimates in Euclidean space[21IS] and on non-trapping manifolds (19) . 

Our method is to introduce an analogue for the wave equation of the Heisenberg equation for the 
Schrodinger equation. There is a self-adjoint operator H which determines the time evolution of the solution. 
We refer to this operator as the Hamiltonian. Other self-adjoint operators are referred to as observables, and 
the analogue of the Heisenberg identity relates the time derivative of a particular inner product involving 
the observable to the expectation value of the commutator between the Hamiltonian and the observable. For 
the Schrodinger equation, the expectation value of the observable is differentiated in time, but for the wave 
equation, the analogue is a more complicated inner product involving both a solution and its time derivative. 
For an operator to be a propagation observable, we require a few other conditions. 

Our first propagation observable is a radial derivative operator directed away from the geodesic surface. 
This is used to prove a smoothed Morawetz estimate, which is like theorem 18.201 but with zero powers of 
the angular derivative and a non compactly supported, but decaying, localisation function. A function, 
of the radial variable, denoted p,, is used to direct this operator away from the geodesies. 

Following this, we define the positive operator L by 

i2 = i_As2, (1.8) 

where A52 is the Laplace- Beltrami operator on the sphere, which acts as the derivative in the angular 
directions. By rescaling the argument of g by powers of L, we generate a new propagation observable, which 
we use to prove theorem 18 . 201 with £ — \- We call this rescaling by L angular modulation. 

To prove the result for all positive e, we use phase space analysis. We introduce the phase space variables, 
which are defined in terms of the original radial co-ordinate, p*, and radial derivative, 

x„ =L"V*, (1.9) 

(1.10) 

op* 
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We localise a propagation observable in phase space by multiplying by a compactly supported function of the 
phase space variable, X{xm) and We also localise with functions which are not compactly supported 

but which decay away from a certain region. By both rescaling the argument of g and localising in 
we can control powers of L in various regions of phase space. However, these estimates require control of 
different powers of L in other regions of phase space. In a process of phase space induction, we are able to 
combine the estimates across all regions, and use the estimate in one region to control the remainder terms 
in another. This is sufficient to prove theorems lOH and 

1.1 Structure of the paper 

Section |21 introduces the wave equation on the Reissner-Nordstr0m solution, some common notation and 
transformations to simplify the equation, and discusses important regions of the geometry, including the 
photon sphere. Section |31 introduces the energy, conformal charge, and a Sobolev estimate. In that section, 
we prove pointwise in time estimates can be reduced to bounding weighted space-time integrals of solutions 
and their angular derivatives. Section ^ translates some of our conditions into terminology commonly used 
by Physicists. In particular, we show that the finite L^, energy, and conformal charge conditions, require 
that initial date, but not its derivatives, vanish at the bifurcation sphere. In fact, functions which grow like 
an arbitrarily small power of the distance from the bifurcation sphere are permitted. 

The remainder of the paper deals with proving weighted space- time integral bounds. Section |S1 provides 
the basic propagation observable frame work, using a commutator formalism analogous to the Heisenberg 
equation from quantum mechanics. Section introduces a propagation observable 7 which majorates a 
decaying weight. This is analogous to a smoothed Morawetz estimate. This is sufficient to control the 
weighted space- time integral needed in the proof of the estimates. The section concludes with an 
estimate for radial data and with a similar result for non-radial data, which reduces the problem to one of 
controlling space- time integrals of one angular derivative in . This localisation is in an arbitrarily small 
neighbourhood of the photon sphere. Section [7| uses a new technique of angular modulation to control | 
angular derivatives in in the desired region. The method is to rescale the previous propagation observable 
by fractional powers of the angular derivative operator. 

In the final section, section |S1 we introduce a family of propagation observables to control 1 — e angular 
derivatives. The propagation observables are like the one from the angular modulation argument, but are also 
localised in both the radial variable and radial derivative. We refer to this as phase space analysis. Subsection 
18. II introduces these observables and rescaled versions of the radial variable and radial derivative, which we 
call the phase space variables. Subsection 18 . 21 provides a commutator theorem which allows us to rearrange 
localisations in the non commuting phase space variables. It also provides lemmas which cover particularly 
common cases. Subsection l8.3l begins the argument by using the commutators to bound fractional powers of 
the angular derivative. These estimates have decaying weights in the rescaled radial variable and derivative. 
Subsection 18.41 removes the weights in the rescaled radial variable, leaving estimates localised in the radial 
derivatives only. All these estimates involve remainder terms which include lower powers of the angular 
derivative, but not localised in the same region. Subsection 18 . 51 combines the results inductively to control 

1 — e angular derivatives with out phase space localisation. It concludes with an L'' estimate with e loss of 
angular derivatives in the energy. 

Through out this paper, the notation C is used to denote constants which may vary from equation 
to equation. Indices are used to separate different constants within a line of an equation and to refer to 
constants later in a proof. 

2 The wave equation on the Reissner-Nordstr0m space 
2.1 The Reissner-Nordstr0m solution 

The Reissner-Nordstr0m solution to the Einstein equations represents the space-time outside a spherically 
symmetric, charged, massive body. It is the unique spherically symmetric, static^, asymptotically flat solution 
to the vacuum Einstein- Maxwell equations, which govern the structure of space-time in the presence of 

^The static condition is redundant since spherically symmetric, vacuum solutions are necessarily static. 
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gravity and electro- magnetic fields. This solution can be represented in terms of the co-ordinates {t, r, 0, cj)) 
by the Lorentzian pseudo-metric: 

ds^ ^Fdt^ - F^^dr'^ - r^idO^ + sin2(6')d(/)^), (2.1) 

where M > is the central mass and Q G M is the central charge, both as measured by observers at infinity. 
The polynomial r^F has two roots at 



r± = Af ± yjAP-Q"^. (2.3) 

In the sub critical case, < < M, these roots determine the location of two event horizons which permit 
geodesies and energy to cross only in the inward direction, towards decreasing r, and not in the outward 
direction. In the critical case, \Q\ = M, there is a single event horizon, which has the same property. In the 
super critical case, |Q| > M, the roots are complex, there are no event horizons to prevent geodesies which 
start at the singularity r = from being extended to r = c», and the central singularity is called a naked 
singularity. We will not study the super critical case. In the non super critical cases, the effects of the event 
horizon and the structure of the interior region are complicated and have been studied extensively [151 \'6'6\ . 
The interior region is called a black hole since anything, including a light ray, may fall in but can not escape. 
We will restrict our attention to the exterior region of the non super critical solution. 



i gK, r > r+ =M + y^M^ - Q"^, {9, 0) eS^. (2.4) 

The famous Schwarzschild solution corresponds to the chargeless case, Q = 0, 

ds^ ^Fdf - F-^dr^ - r^{de^ + sin^ {0)d(j)'^), (2.5) 

r 

t eR, r > r+ =2M, (6*, (f>) eS^ (2.7) 

The Schwarzschild solution, (5 = 0, has only one event horizon at r = 2M, because the second root, r_ 
coincides with the central singularity r = 0. In the exterior region, the Schwarzschild solution is representative 
of all the sub critical solutions. The critical solution is similar; however, F vanishes quadratically instead of 
linearly towards the horizon, and this affects the rate of decay of other quantities. 

In studying the exterior region of the Reissner-Nordstr0m or Schwarzschild solutions it is common to 
introduce a new radial co-ordinate, the Regge- Wheeler tortoise co-ordinate, r, , defined by 

^-F. (2.8) 

This co-ordinate extends from — oo to oo and has the effect of "pushing the horizon to negative infinity" . 
The original radial co-ordinate, r, is now treated as a function of r*. The exterior region of the Reissner- 
Nordstr0m solution is now represented by 

ds^ ^Fdt^ - Fdrl - r^{d6^ + sin2(6')d</)2), (2.9) 

^-1-^ + ^. ("0) 

teR, eM, (61,0) €5^ (2.11) 

In the Schwarzschild case, can be expressed simply in terms of r and M and has simple asymptotic 
behaviour. 

n=r + 2Mlog( — -_)-kC*, (2.12) 

lim r^^r =1, 

hm rrMog f =-i-. (2.14) 



M 

=1, (2.13) 



5 



In the literature, C, is commonly taken, for simplicity, to be — 2Aflog(2) [231121 or — 2Mlog(M/2) ^3]. In 
equation (|2.32() ,we will use a particular choice of C, , coming from the geometry of the Reissner-Nordstr0m 
solution. 

In the general sub critical case , |Q| < A/, the expression for is slightly more complicated, but the 
asymptotic are the same. 



„2 



r*=r + ^ln(— -i) _ln(— + (2.15) 

7'_|- — r_ IVl r+ — M 

lim r^V =1, (2.16) 



lim r:^\nC—^)=-±^. (2.17) 

In the critical case, \Q\ ~ M, the expression for r* and the asymptotics are inverse linear, instead of 
logarithmic, towards the event horizon. 

r — M Af ^ 

=r + 2Mln(— — ) 77 + C'., (2.18) 

Al r — M 

lim rjV=l, (2.19) 

lim r4r~M)=-M^. (2.20) 

r* — ^ — oo 

2.2 The wave equation 

We wish to study the linear wave equation, 

□RN^i=0, w(l)=wo, ^m)=ui, (2.21) 

ot 

in the exterior region of the Reissner-Nordstr0m solution. In terms of the tortoise co-ordinate the d' Alembertian 

is 



□rn -F-\^ - - (2.22) 



(2.23) 



where A52 is the Laplace-Beltrami operator on the sphere. The substitution 

u =ru 

simplifies the wave equation to 

—u + Hu=0, u{l)^uo, —u{l)=ui, (2.24) 

where the operator H is composed of the following terms 

3 

H^J2^-' (2.25) 

1=1 

H.=-^^, (2.26) 
1 dF 

H2 =V{r) = -F—, (2.27) 
r ar 

H3^VL{r){-As2), (2.28) 

Vl(0 =^F. (2.29) 
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We refer to V as the potential and Vl as the angular potential. In the exterior region, r > r+, the angular 
potential has a single critical point at (in terms of the r variable) ^ 



a = ^^^^^^"^ (2.30) 

The critical point, a, is extremely important. Geodesies ai r — a and tangent to this surface will remain on 
the surface forever, and geodesies which approach this surface almost tangentially can approach the surface 
as i ^ oo or orbit the black hole arbitrarily many times before escaping to ±00 |51l32|. This geodesic 
surface at r = a is sometimes called the photon sphere. The region near r = a will also be the most difficult 
in which to prove decay of solution to the wave equation. This region only presents a problem for non- radial 
data. 

The value of r* corresponding to r = a will be denoted a*, and we will introduce the new radial co- 
ordinate 

= r» — a, (2.31) 

This corresponds to taking = with the integration constant C* in equation H2.12|l . (|2.15(l . or (|2.18|) 
chosen so that 



0. (2.32) 



For this reason, = etc. 

For simplicity, we typically use the measure 



(P fj, = dr^d^ fj,g2 . (2.33) 

The space 071 refers to M x 5^ with the measure d'^/i. This defines L'^{Tl) and, more generally, LP{TI) for 
any P > 1- Unless otherwise specified all norms and inner products are with respect to i^(07l). 

The function space iJ^(9Jl) is the collection of functions for which u' and V^^V^u are in i^(9Jl). This 
space is not particularly useful, and we typically consider functions with finite energy or conformal charge, 
as defined in section 13 

When dealing with the original function, u, we also use the measure 

d^fl^r'^dr^d'^Hs^, (2.34) 

and use L^{VJl) to refer to M x 5^ with this measure. 
The two spaces coincide, in the sense that 

\Mmwi) ^\\u\\L^m) (2-35) 
The Schwartz space, S, refers to functions which are infinitely differentiable and for which, for a given 

w e S, 

o j 

Wi,j,ke Z+,3C,,,-fc|p:— V|.i>| < (2.36) 

The Fourier transform on R x 5*^ is defined by first making a spherical harmonic decomposition in the 
angular variable, and then applying the one dimensional Fourier transform on each spherical harmonic. 



3 Methods for pointwise decay 

In Minkowski space, M'^+^, there is a conformal charge which is conserved and which dominates t'^\\r^^W s^uW^ . 
Using a Sobolev estimate, the decay of the angular component of the H^{M.^) norm implies decay of the 
the L^(R'^) norm ^S]. We introduce an analogous conformal charge for the Reissner-Nordstr0m solution. 
The growth of this charge is dominated by the sum of localised space-time integrals of u and its angular 
derivative. Using a Sobolev estimate and the conformal charge, we reduce the proof of a pointwise in time 
LP{VJl) estimate to bounding localised space-time integrals of a solution and its angular derivative. 

^For the critical Reissner-Nordstr0m solution, the angular potential has a second critical point at the event horizon r = M. 
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3.1 Densities, energy conservation, and the conformal charge 

There are many formalisms for studying wave equations. In this subsection, we will introduce various 
densities and show that they satisfy differential relations. Integrating these relations will show that the 
energy is conserved and give an identity for the time derivative of the conformal charge. 

Definition 3.1. Given a pair of functions, {v,w) € S x S, the energy, radial momentum, and angular 
momentum densities are defined respectively by 

e[v,w] =^{w'^ +v''^ + Vv'^ + VlVs'^v ■Vs'^v), (3.1) 

Pp,[v,w] =wv', (3.2) 

Pu,[v,w] =wS/s2V. (3.3) 

The energy is defined to be 

E[v,w] = j e[v,w\d^n (3.4) 

These densities satisfy the following differential relations. 
Lemma 3.2. If u gS is a solution to the wave equation, il + Hu = 0, then the following relations hold 

=^^e[u{t),u{t)] - dp,pp,[uit).Ht)] - Vs2 • (yLpL[u(i),M(i)]), (3.5) 
d 1 

=-g-^PpAu{t),u{t)] - dp,-{u{t)u{t)+u'{t)u'{t) - u{t)Vu{t) - Vs^u{t) ■ VL^sMt)) - • {u'{t)VLV s^u{t)) 

(3.6) 

- ^u{t)V'u{t) - ivs2« • VlVs2U, (3.7) 
and the energy is conserved 



j^E[u{t),u{t)]=0. (3.8) 

Proof. The relations are proven using the method of multipliers, in which both sides of the wave equation 
arc multiplied by a quantity, typically a differential operator acting on u, and then the right hand side is 
rearranged. 

The multiplier gives the relation for the derivative of the energy momentum. 

=(^«)(i2 - u" + Vu- Vl{As2)u) (3.9) 
ot 

=uil + ii'u' + dp,{uu') +uVu + \/s2U- (Vl"^ s^u) + Vs2 • (mVlV52u) (3.10) 
1 d 

=- — {uu + u'u' + uVu + Vs2u • V£,V52u) — dp,{uu') — Vs2 • [uVlV s^u) (3-11) 
d 

=—e[u,u] - dp,pp,[u,u] - Vs2 • {VlPuj[u,u]) (3.12) 

Since the integral of a pure spatial derivative is identically zero, integrating this result gives that 

-^E[u{t), (t)] is zero, and that the energy is conserved. 

The relation for the time derivative of the radial momentum comes from multiplier dp_,u. 

={dp,u){u - u" + Vu + Vl(-As2)u) (3.13) 

=u'ii - u'u" + u'Vu + M'Vr,(-A52)w (3.14) 
d 

=—{u'u) - it'll - u'u" + u'Vu + Vs2ii' • VlVs2m - Vs2 • (■u'VlVs2u) (3.15) 
d 1 

=-g^Pp, — 9p,-{uu + u'u' — uVu — Vs2u • VlV s^u) (3.16) 

- ^uV'u - ivs2« • Viys-u (3.17) 

- Vs2 • {u'VlVs^u) (3.18) 
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□ 

We now define the conformal charge in terms of the energy and momentum densities. In Minkowski space, 
the conformal multipher is found by conjugating the time derivative with the discrete inversion symmetry 
of Minkowski space. The conformal multiplier is then used to define the conformal charge and its density 
by the same process which defines the energy and energy density from the time derivative. The Reissner- 
Nordstr0m solution does not have a discrete inversion symmetry, so we define our conformal multiplier by 
formally taking the Minkowski conformal multiplier and replacing the Minkowski radial variable by , the 
Reissner-Nordstr0m radial variable. 

Definition 3.3. The conformal multiplier, the conformal charge density for a pair of functions (Vjw) e SxS, 
and the conformal charge for the same pair are defined respectively by 

C={f + pl)^^+2tp,dp,, (3.19) 
ec[v,w] ={t^ + pl)e[v,w]+2tp^pp,[v,w], (3.20) 
Ec[v,w] = j ecb,w]d3/x. (3.21) 

The conformal charge density can be rewritten as a manifestly positive quantity. This form is more useful 
for making estimates. 

Lemma 3.4. For any pair {v,w), 

ec[v,w] =^{t- p^f{w-v'f + ^{t + p^,f{:w + v'f (3.22) 
+ \if + pDVv^ + lit" + pI)Vl{Vs^v ■ Ws^v) (3.23) 
Proof. Only the time and radial derivative terms need to be rearranged. 

{t - p^f{w - v'f =t^w^ - 2t'^wv' + t^v''^ (3.24) 
- 2tp^w'^ + Atp^wv' - 2tpW^ (3.25) 
+ pV-2p^W + pV' (3.26) 

Prom this, 

{t - p,f{w - v'f + {t + p^f{w + v'f =2{t^ + pI){w^ + v'^) + Stp^wv' (3.27) 

□ 

In Minkowski space, since the time derivative is the generator of the time translation symmetry, the 
conformal multiplier is a composition of symmetries, and hence a symmetry itself. From Noether's theorem, 
the conformal charge it generates is conserved. Our conformal multiplier is not constructed from symmetries 
and does not generate a conserved quantity. Heuristically, the change of the Reissner-Nordstr0m conformal 
charge should only involve the potentials, since, formally, the conformal multiplier is the same as in M'^+^, 
and the wave equation differs only by the presence of potentials. The potentials appear in expressions of the 
form 2V + pt,V' . We call them the trapping terms. In R^, the analogue of Vl is r~^, and the corresponding 
trapping term, 2Vl + rVj^, vanishes. 

Lemma 3.5. // it G § zs a solution to the wave equation, u + Hu = 0, then 

j^Ec[u{t),u{t)] = j t{2V + p,V')u^d^ii + j t{2VL+ p.Vi){Vs^u-Vs^u)d^ii (3.28) 
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Proof. We multiply the wave equation by the conformal multiplier, Cu, and then apply the relations from 



lemma 13.21 

={Cu){u - u" + Vu + Vl{-As2)u) (3.29) 

= (t^ + pI)u{u - u" + Vu + Vl{-As--)u) (3.30) 

+ 2tp^u'{u - u" + Vu + yL(-As2)u) (3.31) 

= {t^ + pI) [j^e - dp,pp, ~ Vs2 • {VlPu^)^ (3.32) 

+ 2ip* {^PP' ~ ^P*2^^^ + ~ ""^^ ~ Vs2u • VlVs2u) — Vs2 • (m'1/lVs2u)^ (3.33) 

- 2tp, Q^^'" + l^s^u ■ Vl\/s2u^ (3.34) 

Integrating these terms and then integrating by parts in the angular derivatives eliminates the angular 
gradients. 

= y (t' + pD (^^e - dp,pp, - V52 • {Vlp^)^ d'p (3.35) 

+ j 2tp^ [§iPp' ~ ^("" + ~ ^ '^s"'" ■ ^i^s^") - ^s^ ■ (m'VlV52u) j (3.36) 

- J 2tp^ Q"^'" + ^'^s^u ■ VlS/s^uj Sp (3.37) 
= /(i' + P*) (^e - ap.pp.) (3.38) 

+ j 2tp^(^^pp, - dp__,]^{uu + u'u' - uVu-V s^u-VlV s2u)^ (fp (3.39) 

- y" 2tp, Qw^'" + ^Vs2M • T^lV52u^ (3.40) 
This is further simplified by integrating by parts in the radial variable and isolating pure time derivatives. 

= y ^ {{t" + pD + 2tp,pp,) (3.41) 

+ j -2te-2p^pp,Sp (3.42) 

+ j +2p*pp. + 2t]^{uu + u'u' - uVu - Vs2U ■ VLVs2u)d^p (3.43) 

- J 2tp^ Q"^'" + ^^S'" • K'^S^u^ d^f^ (3.44) 
d f 

^Jt I 2tp,ppJ^p (3.45) 

- y" 2t(yu2 + Vs2u • l/LVs2'u)dV (3.46) 

- J 2tp^ Q"^'" + \^S''u ■ Vl\/s2u^ (fp (3.47) 

□ 

3.2 Sobolev estimates 



Our goal in this section is to bound the norm by the energy, the conformal charge, and a negative 
power of t. The main step in this is to prove a Sobolev estimate, which, roughly speaking, controls the 
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norm by the norm. The analogous result in Minkowski space, R'^"'"^, is that the spatial norm can 
be controlled by a third of a factor of the radial component of the norm and two thirds of the angular 
component of the inhomogeneous norm. The energy controls the radial component of the norm, 
and the conformal charge controls the product of a positive power of t and the angular component of the 

norm. It is not obvious from the definition that the conformal charge controls the weighted norm 
needed in the inhomogeneous part of the angular norm, but it does 18 . A similar result holds for the 
Reissner-Nordstr0m solution, but, because of the weight appearing in the norm, we require two estimates 
on weighted norms instead of one. 

The two estimates in the following lemma are the estimates needed for the Sobolev estimate. Note that 
the terms appearing on the right are independent of w, the second argument of Eq. 

Lemma 3.6. There is a constant such that, for all (v, G S x §, 



Ec[v,w]>C{v,^^v), (3.48) 



Ec[v,w]t-^ >C{v,^^v). (3.49) 
pi + 1 



Pi 

Proof. This is a smoothed version of the argument used in M" |18) . 

The ingoing and outgoing wave terms can be isolated and rearranged. We introduce the notation ec,(t,p,) 
to denote twice these terms. 

ec,{t,p,)[v,w] ^\\p,wf + \\twf + Wp^vf + \\tvf + {wAp*tv) (3.50) 
= \\tw + p^v'W^ + Wp^w + tv'W'^ (3.51) 

The weighted term hu = h(p^)u is introduced into these terms. 

^C,{t,p,)[v,w] —\\tw + p^v' — + Wp^w + tv' + thv\\'^ (3.52) 

+ 2{p^v' , p^hv) — 2{tv' , thv) — {p^hv, p^,hv) — {thv, thv) (3.53) 

Dropping the first two terms, which are strictly positive, and integrating by parts in the following pair yields: 

ec,(t,p.) [«, w]>- {v, {2p,h + {pi - t^)h')v) - (v, {pi + t^)h^v) (3.54) 

A particular choice of h can now be made in terms of parameters a > and e > 0. 

/i(p*) —e—^ — (3.55) 
Pi + a 

h'{p.) ^e^^ (3.56) 
We now substitute this choice of h into the previous calculations. 



'pi + a^^P* \pl + aY' {pl+af 



{2p.h + {pi - e)h') - {pi + t^)h' = - e( + [pl t') , , :\ ) - (3-57) 



pt + 3apl - t^a + i^pl 2 Pt + t^pl .q^o^ 
" ' {Pl + af ' {pl+af ^'-^'^ 

I ^ 2n_P^ 3ap^ -t^a 2n t^Pl 



{pl+a)'^ ■'(p2 + a)2 '-(p2_^Q,)2 ^"'"'{pl + aY 

(3.59) 

For e S (—1,0) the first, second, and forth terms are positive. For e S (—^,0) and pl > |, the forth term 
dominates the third by a factor of 2. e can be chosen sufficiently close to zero so that for pl < j, 

{pi + ay {pi + ay 
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Thus, with a = 1, 

Ec[v.w\ > ec,it.p.)[v,w] + {v,H2v) >{-e-€^){v, ^^^^v) (3.61) 

p* + 1 

This proves the first statement in the theorem. The second part follows trivially by dropping the pl/{pi + 1) 
term from the first estimate and dividing hy t^. □ 

We now turn to the Sobolev estimate. Because our main interest is making estimates in terms of the 
energy and conformal estimate, the estimate is expressed both in terms of norms of the derivatives of the 
function under consideration and in terms of the energy and conformal charge. Once again, the energy and 
conformal charge take a second argument w which does not appear in the quantity estimated. 

Lemma 3.7. There is a constant C , such that if {v, ui) G S x S, then 

||F5r^w||i6(OT) <C{E[v,w]+Ec[v,w]t-^)^Ec[v,w]h'i. (3.62) 

Proof. We begin by proving a Sobolev estimate which controls the norm by weighted norms. 

Following the standard argument PHIEIj the proof starts with a W^'^ ^ estimate. Take ^p{p^,uj) 
infinitely differentiable and of compact support for simplicity. To this function ip associate a function on the 
sphere 

f d 

Ii{uj) = J |— V'(p*, w)|dr* (3.63) 

Integration and the Sobolev estimates for M and S"^ can be applied. The Sobolev estimate in 1 dimension 
follows directly from the fundamental theorem of Calculus. 

mp*,u;)\<h{co) (3.64) 
mp*,u;)\^ <hico)^mp.,Lu)\ (3.65) 

\^p{p^,uj)\iduj < Ii{Lu)^\ip{p^,uj)\du; (3.66) 

The second term in the final product on the right side is estimated by the spherical Sobolev estimate. 
The Sobolev estimate on 5^ follows from using a partition of unity on into co-ordinate charts and then 
applying the Sobolev estimate on j^. The notation || V"!! (i) (p*) is introduced to denote the L^{S^) 

norm of 'ijj{p^,uj) with fixed, and the norm of {TS'^ valued) vector is equal to the norm of the length 

\\vs^nx = \\{\vs-m\x. 

^ \i:{p,,uj)\idij <{ [ Ii{u;)duj)^[ \^{p.,uj)fdcj)i (3.67) 
<II^V'l|f (llVs^VI!(i)(p*) + ll^ll(i)(p*)) (3.68) 

A radial weight /"(p*) can be introduced before integrating with respect to dr^. For this proof, the exponents 
a and P are used. They have no relation to the value of r — a governing the location of the photon sphere. 

/ / r(pOIV'(p*,c^)l^dt^dr, <||^V'llf /r(llVs2V'll(iUp*) + ll^llfM(P*))dr* (3.69) 

ii/^V'iii <ii^V'ii!(iirv52V'iii + iirv-iii)^ (3.70) 
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The substitution ip = f^\vi\'^ transforms this to a Hi ^ Sobolev estimate. 

<(4 J f\vA'\^^vi\+Pf-'f'\vi\U'f,f^ (3.71) 

X r+^\vl\^\Vs2Vl\d^^i + J r+^\vifd^f,)^ (3.72) 

X ( j /2"+2'3-l|t;i|6d3^)i( j /|(-As.)^«i|2dV + 1 /IflPdV)^ (3.74) 

We now choose 

a + ^/3=2/3 = 2a + 2/3-l =^ a =i, /? =1, (3.75) 



so that 11/ 3 Will 6 can be cancelled. 



5 f\, „2, 



ll/^wille <4(||^wi||2 + ||ywi||^)3 (11/2 (_A5,)3„,||2 + ||/2«i||2)3 (3.76) 

The weight / = r^^ is now taken to give a result analogous to the Sobolev estimate in M'^. 

IjrT't'ille <C(||^t;i||2 + || - 2Fr~\if)'^\\r-\-^s-)'^v42 + Ik^'^^ilU)^ (3.77) 
op* 

To complete the proof of the Sobolev estimate for the Reissner-Nordstr0m solution, the substitution 
vi — F^v and the inequality Fr^^ < (1 + P*)^^ a-re used. 

||i^^r^w||6 <C{\\F^^u\\2 + \\lF^2Mr-^ - 2Q^)v\\ + \\ - 2Fh-'uf)^ (3.78) 
op* 2 

X {\\Fh-\-As2)K\\2 + llF^r^^wlla)^ (3.79) 

<Ci\\F-2 J^v\\2 + + pl)-h\\)-^ (3.80) 
op* 

X (||F^r-i(-As2)3t,||2 + ||F5r-iw||2)3 (3.81) 

We introduce the dummy function w to act as the second argument of the energy and conformal charge. 
From the definition of the energy, the conformal charge, and lemma 

||^^^r-i(-A52)^«||2 <Ec[v,w]t-^, (3.83) 

\\Fir-^v\\l< + ply^W^ <CEc[v,w]t-^. (3.84) 

Substituting these into equation (|3.81|l proves the result. □ 

We remark that if we repeat the same argument with / = Fr^^ and vi ~ v, the Sobolev type result, 
analogous to equation (|3.81(l , would be 

llF^r^wlle <C(||— w||2 + ||(-2r-i + 6Mr~^ ~ iQ\-^)v\\)^ (3.85) 
op* 

X (||F5r-i(-A52)5„||2 + ||i^^r-iw||2)i (3.86) 

The new weighted norm is also controlled by lemma for < > 1, although, there is no longer the 
additional factor of t^^. 

\\i-2r-' + 6Mr-^ - AQ\~^)vf < ||(^i±^u)3«|| <Ec[v,w]. (3.87) 
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Combining these results, we have control of a norm which decays less rapidly towards the event horizon but 
at the cost of less time decay on the right hand side. 

J \v\^F^r-^d^n <CiE[v,w] + Ec[v,w])Ec[v,w]h-^, (3.88) 
3.3 Local support of the trapping terms 

We refer to terms of the form 2V + p*Vl, for both the potential and the angular potential, as trapping 
terms. In this section, we show that the trapping terms are positive only in a finite interval of p* values in 
the subcritical case. The functions W and Wi will refer to the positive part of the trapping terms. Through 
the conformal identity and the Sobolev estimates, this reduces the problem of finding point wise, weighted 

estimates to proving local decay estimates of the form J J xlul^t^lJ-dt + J J s^'^l^d'^ f^dt < C. 

One of the factors in the derivative of the potential will require careful attention, both here and in 
subsection 16.21 We introduce it with the notation Pgir). 

Definition 3.8. For Q e [0,M], Pgir) is defined by 

Pq (r) = Mlr^ - 4(Q2 + 2M'^)r^ + IbMQ^r - QQ^ (3.89) 

We now show that the trapping term for the potential is positive only in a bounded set of values by 
computing the limit at ±cxd. 

Lemma 3.9. The derivative of the potential is given by 

V' = -2Fr-'^ Pqir) (3.90) 
For sufficiently large, 2V + Pi,V' is negative. 

There is a compactly supported, positive, bounded function, W, such that W > 2V + ps^V 
Proof. The derivative is computed from the definition of V in eauation l2.27l 

d 

=F—((r-^ - 2Mr-2 + Q^r-^)(2Mr-^ - 2Q'^r-^)) (3.92) 
or 

d 

=F—(2Mr-^ - 2(Q^ + 2Af^)r"'^ + eMQ^r^"" - 2QV"^) (3.93) 
or 

= - 2Fr-^(3Afr^ - A{Q^ + 2M^)r^ + IbMQ^r - QQ^) (3.94) 

From this, 

2V + p.V =2r-'F (^^ - ^Pq(o) (3.95) 



r 



To show that this is negative for sufficiently large values of |p, |, we will multiply by a positive factor and 
then show that the resulting quantity has a negative limit as p* — > ±cx). As p* — s- — oo, the subcritical and 
critical cases must be dealt with separately. 

In both the subcritical and critical cases, for p* ^ oo, ^ ^ 1 and F ^ 1, so 

r^{2V + p^V') ^2M -iM ^ -M (3.98) 
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For — > — oo, the original radial variable has limit r ^ r+ = M + \/ — Q'^, and the limiting value 
of PQir+) is 

Pq{M + y/AP - Q2) ^ _ 2(2Af2 - Q2^(Af2 - Q^) _ uiy'AP - Q'^{M^ - Q^) (3.99) 
In the subcritical case, the term with — dominates, so it is sufficient to look at Pq(r^). 

— f^(2F + p.V) ^Pq{t) (3.100) 



= - 2{2M^ - Q^){]VP - Q2) _ 4:My/AP - Q^{M^ - Q^) (3.101) 
<0 (3.102) 

In the critical case, since Pm{M) = , it is necessary to multiply by a different positive factor before 
taking the limit. For p.^, —00, from the asymptotic behaviour of given in equation H2.2()|l 

:{2V + p.V) ^2M—— + -Pl--PM{r) (3.103) 



2Fr~Ar ' r~M r-Mr^ 

3Af 

=2M + p,(r-M) — (r-2Af) (3.104) 

-^2M-3M (3.105) 
= -M (3.106) 

From these limits, it follows that 2V + p^,V' is negative for sufficiently large values of |p*|. Since V and 
p<,V' are continuous, it follows that 2V + p*T^' can be bounded above by a compactly supported, positive, 
bounded function. □ 

In the subcritical case, a similar calculation of limits as p^ ±oo shows that 2Vl + p*V^ is positive only 
on a bounded set of p* values. In the critical case, the angular trapping term is positive as p* —>■ — oo, and 
we compute the rate at which it vanishes as r ^ r+ = M. 

Lemma 3.10. In the subcritical case, for |p*| sufficiently large, 2Vl + p*Vf^ is negative, and there is a 
compactly supported, positive, hounded function, Wl, such that Wl > 2Vl + P*Vl- 

In the critical case, there is a positive, bounded function, Wl, such that Wl > 2Vl + P*Vl, Wl is 
identically zero for sufficiently large p*, and Wl decays like {r — M)F as p* — oo. 

Proof. As in the previous lemma, 2Vl + P*V[ will be shown to be negative from the fact that when a positive 
factor is applied, the limit as p* — > ±oo is strictly negative or the quantity diverges to negative infinity. For 
this lemma, it is simplest to separate the subcritical and critical cases, when evaluating the limits. 
From the definition of Vl , 

-2F ( 3Af 2g2\ 
d,.yL=—r[^~ — + ^] (3.107) 



so that 



2F-p,F'=^(^r-p4l-^ + ^)) (3.108) 
In the subcritical, as p.^ — s- oo, by the explicit expansion of p* given in eauation l2.15l 

[2V - pX) ^r~{r + log(^) - l°g(^) + - — + ^) (3-109) 

2 9 

+ log(— -±) + -— log(— — ) + 0(l) (3.110) 



r_|_ — r_ M f+~ f- M 

oo (3.111 
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In the critical case, as p, — > cxd, by the expansion in 12. 181 

^ {2V - p.V) =r-{r + 2MlogC—^) - ^ + C)(l - ^ + ^) (3.112) 
21' M r — M r 

= -2Mlog(^l^) + 0(l) (3.113) 

-^-oo (3.114) 

For —Qo, we use both the asymptotics of from subsection 12 . II and the geometry of the angular 

potential from subsection 12.21 The critical points of Vl are the roots of — 3Mr + 2Q^, which we denote 

by 



a± = ^ . (3.115) 

Only one of these, a-|_ = a is in the exterior region r > r+. We now multiply the angular trapping term by 
a positive term. 

5 1 3 

^ {2Vl + p.V[) =— +r^- Mir + 20" (3.116) 

lb — p» — 

3 

= — + (r-a)(r-a_) (3.117) 
-p* 

In the subcritical case, the relevant terms are ordered q;_ < r_|_ < a+, so that in the limit p, ^ — oo, r — > r+, 
and 

■^ — {2Vl + p.Vi) ^ (r+ - a_)(r+ - a+) < (3.118) 
lb — p* 

In the critical case, the relevant terms are not distinct, Q!_ — tj^ = M < 2M = a+, so the previous 
argument docs not hold. In fact, the angular trapping term is positive as r — > r+ = M. Instead, equation 
(|3.117|l can be rearranged as 

(r5)(2VL + p.,Vi) =(r3 - p*(r - a){r - (3.119) 
From the asymptotics of p* in equation (|2.18|) . 

/ r — M\ AP 

- p*(r - a_)(r - a+) ^r^ - (r + 2Mln — — — + C*)(r - M)(r - 2M) (3.120) 

\ AI J r ~ M 

= {r^ + A-f2(r - 2M)) + In {^-^^^ - M)2M(2AI - r) (3.121) 

- {r - M){r + C^){r -2m) (3.122) 

On the right, the last term clearly vanishes like r — M, and the second term is negative as p* — > — oo. The 
first term is a polynomial in r which vanishes at r = Af, so it must vanish at least linearly. Thus the right 
hand side is bounded above by a term which vanishes linearly in r — M. The potential term vanishes at a 
rate which is F times faster. □ 

Together, the results from this section show that a weighted norm is controlled, pointwise in time, by 
weighted space- time integrals of a solution and its derivative. 

Proposition 3.11. In the subcritical case, \Q\ < AI , there are hounded, compactly supported functions W 
and Wl such that if u £ S is a solution to the wave equation, u + Hu = 0, then 

\u\%t,p,,u;)F^r-^d^p <C{E[ua,ui] + Ec[u{t),u{t)]t-^)Ec[u{t),u{t)]h^ (3.123) 

Ec[u{t),u{t)]<Ec[uo,ui]+ J J 2t{W\u\'^ + WL\Vs^u\^)d^fidt (3.124) 

In the critical case, the same result holds with Wl zero for r sufficiently large and vanishing linearly in 
r — AI as r —> r^ — AI . 
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4 Relativistic considerations on the event horizon 



To study the behaviour of waves on the event horizon, it is common to introduce the Eddington-Finkelstein 
co-ordinates '^p OSllgSI 

s_=t-p*eM, (4.1) 
s+=t + p^eR, (4.2) 

S"- = - e (_oo,0), (4.3) 

5-+ =e™- e (0,oo). (4.4) 

In fact, this range of co-ordinates only describes the exterior region of the black hole, , and the subcritical 
Reissner-Nordstr0m solutions extend smoothly to positive S- and negative 5"+ ^. The outer event horizon 
corresponds to both the lines S- = and 5"+ = 0. The surface S*. = is the future event horizon, 7i+, and 

= is the past event horizon, 7Y_. To approach 7i_f,, t must diverge to infinite. The bifurcation sphere 
is the sphere given by {S-,S+) = (0, 0), with the spherical co-ordinates free. 

In this section, we discuss how solutions must decay near the bifurcation sphere if they have finite 
conformal charge. 

4.1 Decay on the bifurcation sphere 

To get u e L^{Tl), we made several transformations and change of variables in sectional To determine the 
physical constraints imposed by finite energy and finite conformal charge, we must return to the original 
function u = r^^u, and use normalised vectors and their duals. 



-1 



dr ^F-dt, (4.5) 

dR^F^dp,, (4.6) 

dT=F^dt, (4.7) 

dR=F^dp^. (4.8) 

From this, the natural measure on (/9*,cj) G M x 5^ is 

dfJ-novmaiiscd ^ F^r'^dp^dfis^- (4.9) 

Since there is no bifurcation sphere in the critical case, we will only consider the subcritical cases in this 
subsection. 

We begin by writing the norm in terms of u and the natural measure. 

Mhm = ll«lli2(OT) = \ufr^dr,du; (4.10) 



m 

\u\'^F^'^dHnovm!^lised- (4-11) 

art 



From the first line, at the bifurcation sphere, the function u must vanish at the boundary, but any power of 
f ~ (r — r+) is sufficient, and even inverse logarithmic decay is sufficient. 

^The co-ordinates (s_,s_|_) are typically denoted {u,v) and the term Eddington-Finkelstein co-ordinates properly refers to 
the mixed co-ordinate systems (u, r, 6, (p) or {v, r, 6, 4>). 

*In the critical case, (S_,S+) = (0,0) is a singular point. Otherwise, the critical solution can be extended to an open set 
containing 5- > 0, S+ > and S- < 0, 5- < 0. 



17 



We now write the energy in terms of u. 

E[u, u]= J \u\^ + + Fr-2| Vs2u|2 + 2MFr~^\u\^d^ (4.12) 

= J + \Fr-^u + u'\ + Fr-^\\/s2u\^ + V\u\^y(fn (4.13) 

= J + l^f + Fr-^l^s-^l^yd^fi (4.14) 

+ J 2Fruu + F^\uf + Vr^\uf(ffi. (4.15) 

The co-ordinates t and are singular at the bifurcation sphere, so we must rewrite the derivatives in terms 
of ut — Otu, uh = dRU, and the four gradient \I ^^u. 

E[u,u]= J {F-^\u\^ + F-^\u'\'^ + r-^\\/s2u\^)Fr^dp^d^Lj (4.16) 

+ / -(P^r + P'M^ + + rF^M'd'fi (4.17) 
J dr dr 

= J (l"Tp + \ur\^ + r-^\Vs2u\'^)Fr^dp^d'^uj (4.18) 

\\'4u\'^F^dpnoTm!,\iscd- (4.19) 



For the energy to be bounded, it is sufficient that, far from the bifurcation sphere, the four gradient of 
u is integrable and that near the bifurcation sphere, 

|V4Up <C(r-r+)-i+% (4.20) 

|V4u| <C(r-r+)-5+^ (4.21) 

Under this condition, 

\drii\=F-^\dp,u\ (4.22) 

(4.23) 

<F"5|V4m| (4.24) 

<C(r-r+)-l+^ (4.25) 

The additional condition that u vanishes at the bifurcation sphere restricts the size of u near the bifurcation 
sphere, 

\u\<C{r-r+Y. (4.26) 
From integrating the length of the radial vector Or, the distance from the bifurcation sphere, s, is given 

by 

s~(r-r+)i (4.27) 

Thus, the sufficient condition on the energy in eauation l4.21l Dermits functions which grow like an arbitrarily 
small power of the distance from the bifurcation sphere. 

The calculations and conditions for the conformal charge are similar. Once again, we begin be rewriting 
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the conformal charge in terms of u. 

Ec[u,u]^ j {pl + l){\u\'^ + \u'\^ + Fr-'^\VsM'^ + '^MFr-'^\u\^)(f^x (4-28) 

= jipl + 1)(|wt|' + \UR? + r-^\VsM'')Fr^dp,cPu (4.29) 

+ j {pl + l){2Fruu' + F^\u\^ + Vr'^\u\'^)cPpi (4.30) 

= jipl + 1)(|wt|' + + r-'|Vs2u|2)i^r2dp,d2^ (4.31) 

+ J -2Frp^\uf(ffi (4.32) 

(p2 + i)(|uj,|2 + |^^|2 + r~^\^g,u\2)Fr^dp,d^Lo (4.33) 



^lul^Fr^rfp^d^^ (4.34) 
r 

= J |V4Up(p^ + l)F^d^„ormalisod + ^ | W P (" ^ )-F^ dAinormalisod (4.35) 

The weights in the conformal charge only differ by factors of which is logarithmic in (r — r-|_) near 
the bifurcation sphere. Thus, the sufficient condition given for the energy in equation 14.211 is sufficient to 
guarantee that the conformal charge is also bounded. 

5 The Heisenberg-type relation 

For the Schrodinger equation from quantum mechanics, 

-i^t + H^P^O, (5.1) 

there is the well known Heisenberg type relation for the time derivative of the expectation value of a self- 
adjoint operator, A, 

(5.2) 

This formulation is central to the standard interpretation of quantum mechanics which associates operators 
to physically observable quantities, and the expectation value to the mean observed value. This formulation 
was also used in the original proof of scattering for the quantum n-body problem |12l 1201 1^ I29| . 
We begin by defining the commutator in the form sense. 

Definition 5.1. If H and A are two self-adjoint operators, and D{A) n D{H) is dense in L^(9Jl), then the 
commutator [H, A] is defined to be the form fl given by 

n{v,w)^ {Hv, Aw) ~ {Av, Hw) . (5.3) 

The Heisenberg-type relation follows from the wave equation and this definition. 

Lemma 5.2. If A is a time-independent, self-adjoint operator, and u G § is a solution to the wave equation, 
u -\- Hu = 0, then 

j^{{u,Au) - {u,Au)) =(u, [H,A]u), (5.4) 
where {u, [H,A]u) is understood to mean the quadratic form [H,A] evaluated on the pair {u,u). 
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Proof. We begin by computing the left hand side. 

^((u, All) - {ii, Au)) ={u, Ail) + {u, Ail) - (ii, Au) - {ii, Ait) (5.5) 

={u, -AHu) + {Hu, Au). (5.6) 

Since A is self-adjoint, 

^{{u,Au) - {u,Au)) = - {Au,Hu) + {Hu,Au). (5.7) 

The right hand side is exactly the commutator. □ 

Our method will be to find bounded propagation observables. A simple example of a propagation ob- 
servable, A, which majorates an operator G, is one for which 

[H,A]=G*G (5.8) 

(where G* represents the adjoint of G). If A is a bounded operator on the energy space, then, from the 
Heisenberg-type relation, 

J \\Gufdt < J ^i{u,Au) - {u,Au))dt < 4\\u\\\\Au\\ < CE[u,u\. (5.9) 

The standard definition of a propagation observable is broader than this, and our definition will be even 
broader. 

Definition 5.3. Given a pair of operators, A andG, with D{A)r\D{H)r\D{G) dense m-L^(9Jt), the operator 
A is a propagation observable which majorates G if there is a pair of bounded, non-zero operators, Xi and 
X2, for which 

Xi+X2^Id, (5.10) 
[H, A] =G*XiG + lower order terms, (5.11) 

where "lower order terms" refers to the sum of operators R for which either 
1. R is a bounded operator and for all u gS which solve the wave equation. 



J {u,Ru)dt 



< C, (5.12) 



2. or R is an operator with domain D{R) c D((G*G) 2 ), and for all u e D{{G*G) 2 ), 

(m, Ru) < C{u, {G*G)^-'u). (5.13) 

If A maps S to §, even if it is not self-adjoint, then, since H maps § to S, for v G S, H{Av) — A{Hv) is 
well-defined on §. Since H is self-adjoint, by a similar calculation, 

^{{u,Au) - {u,Au)) ={u,{HA- AH)u). (5.14) 

Since products of smooth functions, integer powers of the radial derivative, and integer powers of angular 
derivatives map S to §, if A is of this form, we can write 

j^{{u, Au) - {u, Au)) ={u, [H, A]u), (5.15) 

[H,A]=HA-AH. (5.16) 
For anti-self-adjoint operators, since multiplication by i commutes with H, we can define 

{u, [H, A]u) ={u, [H, -iA]iu). (5.17) 
Prom this, for an anti-self-adjoint operator, we have the same Heisenberg-type relation. 
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6 Morawetz Estimates 



The goal of this section is to prove bounds on weighted space-time norms of solutions. We begin by intro- 
ducing a propagation observable, 7, show it majorates a weighted quantities, and conclude with a Gronwall's 
type argument to integrate the Heisenberg-type relation. These estimates are proven using a spherical har- 
monic decomposition. We show these estimates have a uniform nature in the spherical harmonic parameter, 
so that we can recover an estimate for general u. In later sections, the contribution from H2 is controlled by 
the local decay estimate, and we can use a uniform multiplier. 

In R^+^, the radial derivative can be used as a propagation observable to control the time integral of 
|u(t,0)p. This can be thought of as a weighted space- time integral with the 6- function as a weight. 
Essentially, our propagation observable is a smooth version of the radial derivative, which leads to a smooth 
weight. This follows P^H). 

The spherical Laplacian is well known to have discrete spectrum. On each spherical harmonic we will 
use / to refer to the standard spherical harmonic parameter and / to refer to the value of the square root of 
the Laplacian on that harmonic, hence r = l{l + 1). We will also introduce the effective potential on each 
spherical harmonic, given by 

Vi=H2+H:i = V + ?Vi. (6.1) 

If necessary, we will use P; to denote the projection on to the spherical harmonic with parameter I, but 
typically, in this section, we will assume that u has support on only one spherical harmonic. In this case, u 
satisfies 

u^u" + PViu = 0. (6.2) 

We will show in lemma IHl^ that Vi has a single critical point, which is a maximum. We denote this critical 
point (a;)*. 

The Morawetz-type operator is defined on each spherical harmonic in terms of the radial co-ordinate . 
The definition involves a weight, g^, which is defined for > 1 so that (70- remains bounded. 

Definition 6.1. Given a > I, b > 0, and I e N, the Morawetz-type multiplier ^i^a is defined by 

giAM^l (TTmF^" ^^-^^ 

1 d d 

=m,bM{p*) + ]^g'ub,AP*) (6-5) 

In all cases a will not vary so the notation 7 — ji^a- o-nd g = gi^„ will be used. The projection operator onto 
the l*^ spherical harmonic is denoted P/. When working on more than one spherical harmonic, we also use 
7 to denote 

^7;,mP;- (6-6) 

We note that as I varies, the observable ji^b.a is simply translated. We will also show that (a;)* 
from which it follows that gi,b,a has a limit in L°°. 

6.1 Preliminary bounds 

Lemma 6.2. If a > 1 and b > 0, then there are constants Ca and Cb, independent of I, such that for all 

{■U,ll,b,aU) =0 (6.7) 

\hu\\L2 <C^\\u'\\l2 + CbWil + pI)-^u\\l2 (6.8) 



<C^y/E[[;] + Cb\\{l + pl)-^u\\L2 (6.9) 
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Proof. Wc work on a fixed spherical harmonic. Since w G S, derivatives can be moved about freely. 

1 ^ d 

(w, ju) ={u, 2 (ff^ + (6.10) 

=0 (6.12) 

The last equality holds since the inner product is symmetric on real valued functions. 

Since g — /q''^''* ^"'^'-'(l + |r|)^'^dr, g' is bounded by Ch(l + pl)^'^^'^ ■ Since the integrand is continuous 
and decays as r^"' for cr > \, g is bounded by some constant Ca- Using this and the fact that u S §, the 
following holds. 

\\^u\\^\gu' +]^g'u\\ (6.13) 

<||.g«'|| + ||ig'7i|| (6.14) 
<a||u'|| +a||(l + p2)-u|| (6.15) 

The constant Ca is given by the integral of (1 + |r|)^'^ from to oo. This is independent of I. The constant 
Cfc is based on the equivalence of 6(1 + — {ai)^\)^" and (1 + pl)^"^'^ ■ Since (a/)* converges to zero by 
lemma 1^31 this equivalence is uniform in L □ 

6.2 Computation of Morawetz commutators 

We show that the propagation observable, 7, majorates (1 — p'i)^'^^^ ■ We will prove this on each spherical 
harmonic, and then show there is a uniform lower bound on these estimates, to get an estimate for uniform 
I. 

We first show that on each spherical harmonic, the eifective potential has a unique maximum. 

Lemma 6.3. For / > 0, M > and < |Q| < M , the effective potential Vi has a unique critical point, 
which is a maximum. The maxima, (a;)*, converge to 0. 

Proof. The Reissner-Nordstr0m potentials are 

I M 9 M n2 

(6.16) 

(6.17) 
(6.18) 

Our goal is to show that in the outer region, the effective potential Vi has a single critical point. The 
outer region is given by r > 7'+ with 



V : 


2M , 2M 
=—(!- — 


+ 91 


Vl 


1 , 2M 
^^'^^~ — + 


€■ 


Vi 


^V + PVl. 





r± =M ± a/M2 - Q2^ (6.19) 

being the zeroes of the factor 1 — 2M /r + Q^/r^. 

We will do this by showing that a positive factor times the derivative of each of the potentials increases 
from a nonpositive value to infinite from r — 2M to r — s- cx) and that the potentials are negative for r from 
the horizon at r = r+ to r = 2M. 

The derivatives are given by 

V' =- 2Fr-''P, P =3Mr^ - 4:{Q^ + 2M^y + IdMQ^r - 6Q^, (6.20) 

V[=- 2Fr''^PL, Pl =r^ - SMr^ + 2QV^ (6.21) 

V/ =y' + Pvl -2Fr-'^Pi, Pl =P + PPl. (6.22) 
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We will also use 



/ =r-^P = 3Mr - 4(Q2 + 2M^) + 15MQ^r~^ - GQ^r'^, (6.23) 
lL=r-^PL=r^ -iMr + 2Q^, (6.24) 
Ij^r^^Pi^I + FlL. (6.25) 

Showing that Pi has a single root in the exterior region is equivalent to showing that Ij has a single root in 
the exterior region. 

We start by showing that / and II are increasing in the region r > 2M. The derivative of / is 

dri =3M - 15MQ^r-^ + 12Q\-^ (6.26) 
=r"^(3Mr^ - 15A/QV + 12Q'^). (6.27) 

The derivative of the term in brackets, for r > 2M and Q < M, can be estimated by 

driUdr^ - 15MgV + 12Q'*) = 9Mr^ - 15MQ^ > 0. (6.28) 

At r = 2M, the quantity in brackets in I^17\i is 24M'^ - SOAPQ^ + 12(5** > 0, thus drI is positive, and / is 
increasing from r — 2M . The derivative of 7^, driL = 2r — 3A/ > 0, is positive for r > 2M (and even for 
r > 3M/2). 

The value of these quantities &i r — 2M is 

/ =6A^f2 - 4(Q2 + 2M2) + 7.5g2 - (6.29) 

2Af ^ 

^^{-2M^ + 3.5Af2g^ - l.SQ'^) < 0, (6.30) 

Il =4M2 - + 2Q2 < o. (6.31) 

Thus each of / and II are increasing for r > 2M, diverge to infinite as r —^ oo, and are nonpositive for 
r = 2M. 

We now show that / and II is negative for r between r+ and 2M . To show / is negative, it is sufficient 
to show P = r'^I is negative. This is a cubic, and we observe a number of properties about it before drawing 
a conclusion. The derivative is quadratic with two roots, 

dr{P) =9Mr^ - 8(g2 + 2A/2)r + ISMg^, (6.32) 



4g2 + SAP ± Jmm^ - 7iM2g2 + i6g4 

f = 7:t7 • (6.33) 

The discriminant is bounded above by 64M2 - 64Af2Q2 _^ ^qqA ^ ^g^j2 „ 4g2)2 xherefore, the upper root 
is less than 16Af/9 < 2M . The discriminant is bounded from below by 9g'*, so the lower root is bounded 
above by (SM^ + Q'^)/9M < M. The value of P = Pi at r = r+ is 



P{r+) = - 2{2AP - g2 + 2M^/IP~~Q^){AP - Q^) < 0. (6.34) 

Considered as a function on the entire real line, P goes from negative infinite to a local maximum value at 
the lower root of the derivative, and then decreases to a non positive value at r = r+. Since Pi is cubic, 
it can only have one root for r > r^. Since that root is is in r > 2M, the cubic must be negative for 
r e [r+,2A/]. 

The function II is much easier to estimate. It has two roots at 



3A//±v/9Af2-8g2 
r= . (6.35) 

For g < M, the lower root is less than or equal to M and the upper root is greater than or equal to 2M. 
Thus, Il is negative from r+ = M + y/Q'^ - NP to 2M. 
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We know that / and II are negative for r between r+ ~ M + ^ _ jy{2 g^j^^j 2M, are increasing for 
r > 2M, and go to infinite as r — > oo. Therefore, the weighted sum = / + /^Jj;, has the same properties 
and must have a single root for r > r+. 

Since the derivative goes from negative to positive, these critical points are maxima. Since Vil^^ Vl 
and the maximum of Vl occurs at 0, if (a;)* denotes the maximum point of VJ, then (a;)* ^0. □ 

We now show that the commutators uniformly bound a polynomially decaying term, for an appropriate 
choice of b. This value of b will always be used in 7;^^ from now on. There are additional positive terms 
which are also dominated, but we are better able to take advantage of them in the later sections using a 
slightly modified multiplier. 

Lemma 6.4. If a > 1, M > 0, and \Q\ < M, then there is a choice of b for which there is a constant C 
such that for all m G S, 

(u,[iJ,7]) >C{u,{1+pI)-^-'u). (6.36) 

Proof. We work on a single spherical harmonic and sum at the end of the argument. The commutator can 
be computed exactly as 

We will use x to denote — (a;)* 

g' =- — 5—— (6.39) 
-&Vsgn(x) 

^ {i+b\x\r+^ ^ ' 

To estimate the contribution from the S function, we introduce a smooth, weakly decreasing, compactly 
support function, % which is identically one in a neighbourhood of x = (this is unrelated to any other 
cut off function used elsewhere in this paper). Using integration by parts, and temporarily treating u as a 
function of x, 

g^{xxu^)dx (6.42) 

= J xu^dx + J xx'u^dx + J 2xxu-^—udx, (6.43) 

J xu'^dx < J {\xx'\ + x^xWdx + J xi-^ufdx, (6.44) 
f°° d 

u{Of=- —{xu'^)dx (6.45) 
io 



d_ 
dp* 

d_ 
dp. 



f , 1 f 
= - / xu dx- 2xu 

Jo Jo 



udx (6.46) 



roo POO roo q 

< I \x!\u'^dx+ I xu'^dx+ / xii^ufdx. (6.47) 
Jo Jo Jo op* 

Applying a symmetric argument for (—00,0] and substituting the result for _/ x'^^dx, we have 

2u{0f < J (Ix'l + \xx'\ + x\)v?dx + j 2x{^ufdx. (6.48) 
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We use a similar integration by parts argument to estimate (l + 6|a;|) " 

d 



= y -^{g"v?)dx (6.49) 
= y g"'u^dx + J 2g"u^udx (6.50) 
'^^'^^^^ rrfa;<2a6V0)' (6.51) 



(1 + 61x1) 



cr+2 



^ 7 — ,, ,^ ' M^rfx / b , "7, dx (6.52) 

(l + 6|a;|)'^+2 y W fj+l (l + 6|a;|)'^ i ^ ^ 



Substituting equation 

(1 + 6|x|)^'-+2 ^^ -^^'/ + \xx'\+x\Wdx + ab' J "^Xi^ufdx 



(6.53) 



+ 1 (b^ "(" + ^) + l f^b^^^^dx (6 54) 

+ 27 (l + feN)-+2" '''^+2 7 a + l^l + 6|x|)- ^"-^^^ 

'^^'^ + <2(t62 f (i^'i + + x\)u^dx + 2ab f (^—7^ ItTT + ^bx)i^ufdx (6.55) 



{l + b\x\y+^ - 7''-^' J '(T + l{l + b\x\y '^dp^ 

The commutator term can now be estimated 



'"u^ - gVi'u^dx (6.56) 



2^ -{4-ufdx + a&2u(0)2 _ 1 / _^M^±il_u2rfa; - gVlu'dx (6.57) 



(l + 6|x|)'^'ap» ' 27 (l + 6|a;|) 

^/(TTW^^")'^^ + ^^'"(°)' (^-^^^ 

- [ (Ix'l + Ixx'l + x\)u'dx -<jb [ { l^,^ J., + 2bxK^ufdx - gV/u'dx (6.59) 
J J a + I [I + b\x\)'^ dp* 

>2b J - ^^^)(^")'^^ + / (-'^^'(Ix'l + I^X'I + ^'X) - 5^/)w'da: + ab^Of 

(6.60) 

The term (t6m(0)2 is positive and wiU be ignored. The term {a + 1)^^ {1 + b\x\)^'^^'^ — abx is independent of I, 
so we can choose b sufficiently small so that this is bounded below by c(l + b\x\)~"'~^^. Since V/ = V + V^P, 
and the (a;), converge, the value of V/' is uniformly bounded from below at (a;)*, and b can be chosen 
sufficiently small so that the quadratically vanishing quantity —o'b'^{\x'\ + \xx'\ + x^x) ~ 9%') is bounded 
below by —cx'^X: uniformly in I. Thus there are uniform constants, ci and C2, for which, 

J 25'(^")' - l^u' - gVydx (6.61) 

^^2/ JT^^^^'dx. (6.63) 

Since the (a;)* converge, there is a uniform equivalence between {l+b\x\)~'^~'^ and (l+p^)-'^/^-!. Integrating 
over the angular variables gives 

{u, [H, -i]u) >C{u, ——^—^u). (6.64) 

(1 +p2)2+^ 

□ 
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6.3 local decay estimate 

We begin with a type of GronwaU's estimate. 

Lemma 6.5. 1/9 : [0,cx)) — > [0,cx)) with uni/ormly bounded derivative, e £ (0,^), and there are constants 
Ci, 6*2, and T , such that for t > T 

f 9{TfdT < Ci + C2fe{tf~' (6.65) 

then there is a sequence {ti} such that Unii^oo ti — oo and 

hm tlOitf-' = (6.66) 

i — *CXD 

Proof Successively stronger bounds on 9{t) will be proven. 

The bound on the derivative implies 9{t) and t'^9(t)^~'^ are linearly bounded above. 

Suppose there is no sequence {ti} on which 9{ti) — > 0, then 9{t) is bounded below by a constant. By 
the integral condition t'^9{t)^~'^ is bounded below by a linear function. From this 9(t)'^ is bounded below by 
a quadratic, and its integral is bounded below by a cubic. This contradicts the linear upper bound of the 
integral. Therefore, there is a sequence {ti} such that — > oo and 9{ti) 0. 



Since e < ^, there exists r < 1 such that 



1 -e 

< 



2 - r + re 1 - e 
Now choose 5 negative such that S > from which it follows that 

9 A 

r6 < 1 + (6.67) 

1 — e 

Suppose 3K > 0, S* > : Vt > S* : 9{t) > Kt^ , then by the integral condition H6.65|l there is a positive C3 
such that 

6*3*^+2* < Ci + C'it'eity-" (6.68) 

Since S > > 1 + 25 is strictly positive. Thus for sufhciently large t, there are constants C4, C5 such 
that 



C,t'+^ < 9{t) 



If 5 is sufficiently close to zero, then this contradicts the previous result that 9{t) goes to zero on a subse- 
quence. If 5 is larger, then we use the fact that by equation H6.67|) . 9{t) > C^t^^ . This implies the original 
assumed lower bound is replaced by the larger lower bound i'"*. Repeated iterations of this process shows 
that 9{t) is bounded below by t^"^ for any n. Since r < 1, this reduces the situation to the 5 close to zero 
case, which led to a contradiction. Thus 9{t) can not be bounded below by a function of the form Kt^ for 
5 > In particular, 9{t) can not be bounded below by a function of the form Kt^^ and so the desired 
subsequence must exist. □ 

The propagation observable 7 majorates the weig ht {l + pl) ' for /? > |. This is a regularised version 
of the Morawetz estimate. As a consequence, the solution is space-time integrable in any bounded region, 
and must decay in that region. 

Theorem 6.6 (Local Decay Estimate). If (3 > |, then there is a positive constants K^d and such that 
for u £ S a solution to the wave equation, u -\- Hu = 0, 



where E — E[uo,ui]. 



\\{l + pi) — urdt <Kld^E{^/E+\\u{1)\\) (6.69) 
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— p 

Proof. First a crude estimate on the norm growth of || (l + p^) ^ m|| is made. In this paragraph, we use the 
notation f ee [l + pi) ^ < 1. 

^^\\fu\\' ^^^{fuju) (6.70) 

2\\fu\\^^\\fu\\^2(fu,fu) (6.71) 
d 



g-^WM <\\fu\\ < \\u\\ < VeW)] (6.72) 

Initially (3 is restricted to (|, 2), and the variable a — 2(/3 — 1) G (1, 2) is used. This a is distinct from the 
a used previously in the commutator argument and definition of 7. From the commutator estimates lemma 
16.41 there is a constants C2 such that 

^ C2\\il + pl)^ u\\^dt<- {u,[Y^H,,j]u)dt (6.73) 



i=l 

From integrating the Heisenberg relation, the integral on the right is bounded by 

rT Q 

<[^(u,7^^)-2(^,7«)]Li (6.75) 

<[-2(u,7u)]Li (6.76) 
<[||^llll7"ll]*=i + [ll^llll7^^ll]*=T (6.77) 

From energy conservation and the bound on ||7u|| in lemma these norms are bounded by 

[||?i||||7"l|]t=i + [||«||||7"ll]t=T <CiVe{Ve+ II (1 + p2)^u(i)||) (6.78) 

+ C2Ve{Ve+\\{1 + pl)"^ u{T)\\) (6.79) 

Since a S (1, 2), q can be chosen in + 5, f ) C (1, |). If p is the conjugate exponent to q and k = |, 
then 

->l-\ = \ (6.80) 
p 3 3 

(^)9=1 (6.81) 

cr + 2 , ^ 

aq >^ (6.82) 

Holder's inequality can be applied to the last term in equation (|6.79(l using conjugate exponent p and q. 
Continuing to use || • || as the norm gives, 

\\{l + pl)-^ur^f ffl^d^M (6.83) 

<(/ \urd'f,)Hf TTT^rfV)' (6.84) 
Jm Jwt U + P*) 

<HH^H^II '"'T'^ ll' (6.85) 



The fact that aq > can now be used. 



II (1 + Piy^' u\\ <\\u\\H (1 + p^) "V^ ^||i-^ (6.86) 
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Since < E, the norm of u is controlled, < +ti<j2 . The computations from the beginning 

of 16 . 731 to 16 . 861 give an integral inequality for the weighted norm of u. 

T 

II (1 + pD^ufdt < C,iE[u] + Wuf) + C2E[u]^ti\\ (1 + pl)^u\\'-i (6.87) 

Since ^|| (l + p^) " it|| is uniformly bounded by equation (|6.72() . and since ^ G (0,^), the Gronwall's 
estimate lemma 1^751 applies, and there is a subsequence on which 

ti\\{l + pl)^ u\\^-i ^0 (6.88) 

Using the calculations starting from the left hand side of equation H6.73|l , which is monotonically increasing 
in time, to equation H6.86|l . which is sequentially decreasing, the desired result holds for (3 g (|,2). Since 

(l + pj) ^ |u| is monotonically decreasing in /?, the Lebesgue dominated convergence theorem extends the 
result to aU > |. □ 

To conclude this section, we apply this result to the conformal estimate. Since the weights in the local 
decay theorem will dominate any compactly supported function, if u is a radial function, the trapping terms 
can be integrated in time to control the growth of the conformal charge and the weighted norm. 

Corollary 6.7. In the non super critical cases, \Q\ < M, if u ^ S is a solution to the wave equation, 
u + Hu ~ 0, and uq and Ui are radial functions, then there is a constant, C, depending only on ||mo|P, 
E[uq,ui], and Ec[uo,ui] such that, 

\^ 

\u\^{t,p,,uj)F^r'^d^fi)\ <Ct— (6.89) 

Proof. Since the initial conditions are radial, the solution u will remain radial, and 

J T^L|Vs2upd/^ = 0. (6.90) 

Since the potential trapping term, W, is compactly supported, it can be dominated by the weights in the 
local decay result. This provides a bound on the conformal charge. 

Ec[u{T),u{T)]<Ec[uo,ui]+ J J 2TW\u\''d^fidT (6.91) 

<Ec[uo,ui]+t J j2W\u\^d^pdT (6.92) 
<Ec[uo,ui]+tE^E^ + \\uo\\) (6.93) 

From proposition 13 . 1 11 

|M|6(i,p,,w)F3r-'*(iV<C^(^^[«o,"i] +£^c[^i(0,wW]i"')^^c[u(t),M(t)]i^^ (6.94) 
<Ct~2 (6.95) 

□ 

7 Angular Modulation 

To close the conformal estimate it is necessary to estimate one angular derivative in localised near the 
photon sphere r = a. Although this will not be possible, we will be able to bound fractional powers of 
the angular derivative. To do this, we introduce a propagation observable, which is an analogue of 7 but is 
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rescaled, on each spherical harmonic, by a fractional power of 1 — As2. In this section, we control (1 — A52) 2 
in expectation value. 

Inspired by, but deviating from, the standard notation from Quantum Mechanics, L is notationally used 
to refer to the operator square root of 1 — A52; however, since an explicit decomposition of the function 
space into eigenfunctions of 1 — A52 is known, at this point it can still be thought of as a purely notational 
convenience. A strictly positive version of — A52 is used so that inverse powers can be taken in section |H1 
This differs slightly from Z, which was used in the previous section but will not be used again. 

Definition 7.1. L is defined to he the operator square root of 1 — A52. On each spherical harmonic this 
acts as multiplication by y^l + l{l + 1). 

7.1 Angular Modulation and Initial Estimates 

A new operator 7^™ is introduced to give better estimates near ~ 0. Because the contribution from 
H2 = V involves no powers of L, the local decay result allows us to ignore it, and the new multiplier is 
centred at the peak of Vl uniformly in L. Instead, the rescaling parameter will vary with L™. 

Definition 7.2. The angularly modulated multiplier jj^m is defined for m G [0, 1/2] , b > 0, and a > I, by 

gL^ ^ (1 + |r|)-^dr (7.1) 



00 Ml+lH+l))^ p, 

= E/ (l + |r|)-drPz (7.2) 











7.™=i(5L™^ + ^.9.™) (7.3) 

We will take a > 2 since there is no improvement in varying a and it simplifies the estimates. 
Lemma 7.3. If m e [0, 1/2], a >2, and b> 0, then for all u £§, 

IhLr^Ml' < CmiE[u] + II (1 + ply' u\\'). (7.4) 

Proof. Since a > 2, g^m is bounded in absolute value by Co-. The derivative is bounded on each spherical 
harmonic shell. 

g'L^ = L™(1 + feX^lp^l)-" < L'"(l + 6L"|p,|)-2 < C{L"'F-h-^ + (l + piy') (7.5) 

(this constant is independent of the spherical harmonic shell). The norm ||.9^,tiu|P is bounded by ||LF^r^^u|p4 
II (1 + pi) ^ '^'IP which are bounded by the energy and local decay norm. 

WlLr^ur <2{\\gr.^u'\\' + l\\g'^^ur) (7.6) 
<Cm{E[u] + \\{l + piy\r). (7.7) 

□ 



7.2 Direct Angular Momentum Bounds 

The family of multipliers 7^™ is not sufhcient to prove a local decay estimate for the angular energy 
{u,WlL'^u). Instead only an estimate for (u^WlL^u) can be proven directly without phase space anal- 
ysis. The same method as in the proof of local decay is used, and the best result is found by optimising over 
m. We track some of the additional positive terms more carefully to gain this improvement. 

Whenever a commutator involving a localisation or its derivatives, appears in a calculation, that term 
can be controlled by the previous local decay result. For this reason, we start by ignoring the contribution 
from H2 and estimate it later. 
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Lemma 7.4. If m £ [0, a >2, e> 0, and a continuous, compactly supported function x, then there is a 
choice of h, constants Ci , C2 and C3 and a smooth, compactly supported function x such that for all u G S, 

r m 

{u, [H, + i/3,7L"» >Ci{u', -— u') (7.8) 



+ (1 - -){u, -gViiL' - l)u) -C2] b'pixu'dp,, (7.9) 

>e J L^'^xu^dp, + (1 - e) j -qVHL^ - l)u^dp, - C3 j h^plxu^dp,. (7.10) 

Proof. We follow the same argument as before, but with the extra parameters from scaling by L™. From 
the definition of g, we have 

rrriL 

' (7.11) 



(1 + L™&|p*|) 
^(l + L'»fe|p,|)'^+i 



„ -g-L^^fc^sgnp, 



r ^ - 2L-6^.^(,.) + .(. + 1)L-- (,^^,4.|)... - (^-13) 

Once again, we begin by controlling the S function appearing in g'" . For this argument, x denotes a 
compactly supported function, identically one in an interval of width i~™ about the origin and supported 
on an interval of the same scale. 

— {p,XU^)dp, (7.14) 

Xu^dpt: + J pt:x'u^dpt: + J 2p„xu-^udp„ (7.15) 

We isolate the first term, note that the derivative of a characteristic function is bounded by L^™p* times a 
(different) characteristic function (also denoted x), multiply by L'^™, and apply Cauchy-Schwartz, 

L^^'xu'dp, < I L''"'plxu^dp, + 2 (^j L^^ylx'u'dp?! ' (^J L^xi^ufdp,^ ' . (7.16) 

We control the value of the function at p» = using a similar method, integrating only on the positive axis 
to get the value at the origin. 

u{Qf^- ^{xu')dp. (7.17) 

poo Q 

x'u'^dp^ — I 2xu— — udp^ (7.18) 
Jo op* 

<C (^J xu'dp?j ' (^J L^^plxuHp. + j x{^ufdp)j ' . (7.19) 
We now multiply by L to the power 2m. 

L^^u{Qf <C (^j L^"^xu^dp}j ' (^j L^'^plxu^dp, + j L^x{^ufdp)j ' . (7.20) 
Applying equation H7.16|l to control the weighted norms, we have, 

L2'"„(0)2 <Co( j L^^plxu^dp, + j L^x{-^ufdp,). (7.21) 
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We now turn to controlling non compactly supported weights coming from g'" . 

— {g"u^)dp, (7.22) 

= J g"'u^dp^+ J 2g"u^udp^ (7.23) 

"^"^^ / (TTW^^^* <2aL^-bM^f (7.24) 



(i + 6L"'ip*l)'^+2 y i^i CT + 1 (i + 6L'"|p,i)'- y 

(7.25) 

Substituting equation ifT?^ and noting x < C(l + 5L™|p*|)'^, 

/ ir-^T-^Ar^^^dp. (7-26) 

7 (l + 6L™|p,|)'^+2 ^ ^ ^ 

<abCo J bL^-^plxu'dp. + j ^^"^ (rr5^;r^^P* J (7-27) 

, 1 / f^3,»^3 ^(^ + 1) ^J2. , f 2a 2(^u)2 \ 

+ 2 ^ ^ (l + 5L-|p,|)-+2" '^^^ + i ^ (TTfeZn^ ^7 ^ ^ ^ 

/ 7^4^T-^Tr^"''^P* (7-29) 

y (i + 6i'"ip,i)'^+2 ' 

<2abCo J L'^^bplxu^dp. + (^ + 2abCo) J ^^^YT^^^l^'^^* ^^'^"^ 
The commutator term can now be estimated 

25'(^«)' \g"'u^ ~ gVLiL' - iWdp^ (7.31) 

> / ^"fc , ^r, + i2"6Vu(0)2 (7.32) 
7 i + t)|p*| 

- / ^"'"'ranF^"''"" - / - 

> y - 2.5C„)L™6^^^^dp, + L^-b^auiOf (IM) 
aL'^"'b^Coplxu^dp,- [ gV[iL^ ^l)u^dp,. (7.35) 



We choose b sufficiently small so that 2/{a+ 1) — 2abCo > 1/2 and — (e/3)p*V^x > ab'^plx- This choice of b 
can be computed exphcitly from Co and V^(0), but it will be independent of the spherical harmonic. Since 
g vanishes hke L™p* in the support of x, under these conditions, —{€/3)gV[L'^ > aL^^^b"^ pix- This gives us 
sufficient control, except in the case — 1 = when we must subtract off an additional localisation term. 
With these choices, we have 

2g'{^uf - \g"'u^ - gVi{L^ - \)uHp. (7.36) 
op* 2 

^SP' ' Ar. _L ^r2m;,2, ,/n^2 _ /I _ / „T/'rr2 _ 1V,,2^„ _r< / f,2 2, 



>c y ^"^(T^il^Jr^^^* + ci""6'«(0)2 ~ " 3 V ^^^^^ " " ^^Plxdp.. (7.37) 

We ignore the u(0)2 term and divide the remaining positive terms into three pieces. A large piece will be 
1 — e times —gV'^{L? — 1). Another piece will be the (u')^ term and e/3 times —gV'^{L? — 1). In a compact 
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set of width by equation (|7.16|) . this controls L^"^, up to a compactly supported term coming from 

the — 1 — case. The remaining e/3 factor of ~gV[, we use outside \p^\ < In this outer region, 

but inside a compact set independent of L, —gV^{L^ — 1) dominates any continuous, compactly supported 
function, Xi by a factor of > L^™, again with an extra contribution required for the i'^ — 1 = case. 

Here, we use the assumption m < 1/2 again. From this, 

j 'ig'i^uf ~ \g"'u^ - gViuHp^ (7.38) 

>ce y i^^xw^rfp* + (1 - e) y" ffV^i^ - \)v?dp^ J b^plxdp*. (7.39) 

Recall X denotes an arbitrary continuous, compactly supported function, and x is supported on \p^\ < 
CL-". □ 

The next lemma says the contribution from H2 is a local decay term. 

Lemma 7.5. // m e [0, 1/2], a >2, and b > 0, then there is a positive constant C such that for all u lE S, 

\{u, [i?2,7L"»| < C{u, (1 + (7 40) 

Proof. From lemma the derivative of the potential is — V = 2Fr^'^ PQ{r) where Pq is a cubic polynomial 
in r. 

Since for — > 00, r ^ p*, and for p^ — > — cxd, _F decays exponentially, there is a constant C so that 
\V'\ < C (1 + pI) . Since I^l-"! is bounded by a constant and [i?2,7L'"] = — ffi™^', the result holds. □ 

The previous two lemmas shows that is a propagation observable which majorates powers of L. That 
is, ignoring terms which are space- time integrable, the commutator is bounded below by the product of 
powers of L and localisation functions. These localisation functions are functions of L™p^,. In the following 
theorem, the Hcisenberg-type relation with j^m is integrated in time. The theorem shows that is a 
bounded operator on the energy space. 

The localising function x is included to restrict attention to the region bounded but on an L indepen- 
dent scale. In the theorem, taking the optimal value to = ^ gives an estimate for {u^xL'^u) . The estimate 
for other values of to will be used later, in the phase space analysis. 

Theorem 7.6. ^ // to G [0, then there is a constant C such that for all u £ S which satisfy the linear 
wave equation, it + Hu = 0, 

T m, 

/) + {u, ^gV[iL' - l)u)dt <CiE[u] + ||^.(1)||2) (7.41) 



(1 + (L-p*)2 

{u,Lixcu) <C{E[u] + \\u{l)\\') (7.42) 



Proof. Applying the Heisenberg identity to the operator gives 



d ^ 

- ( {u, 7^™ u)-{u, ii))^ {u,Y,H^u) (7.43) 



Using lemma [731 and the fact that 7^™ is antisymmetric with respect to the L'^{dJV) norm when acting on 
C°° functions, this equation can be integrated to get an estimate on the commutators. 



fin, 
•^1 i=i 



jL^]u)dT <2\\u\\\\-fL,^u\\\l (7.44) 

<c{E[u] + II (1 + piy' u{i)r + II (1 + piy' u{t)r) (7.45) 



^This is the equivalent of theorem 7.8 from the previous draft. 
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By the local decay result, theorem 15.61 and lemma [731 

/oo 

II (1 + ^||2^^ < ^(^[^] ^ ||u(l)|P), (7.46) 

/oc 

II (1 + uW^dT < C{E[u] + ||u(l)||2), (7.47) 

(7.48) 



From the first estimate for [Hi + H3, 7l™], equation H7.9|l . the foilowing estimate holds on a sequence of 
times: 



<C{E[u] + \\u{l)r + II (1 + pI)-' u{t)f) (7.50) 

+ CiE[u] + \\{l + pl)-\{l)r) (7.51) 

Since the left hand side is monotonically increasing and the right hand side is decreasing on a sequence, 
by the local decay estimate, theorem (|6.6|) . equation (|7.41|) holds for all t. 

Equation (|7.42() is proven by taking m — ^, and applying the same argument using the estimate in H7.16|l . 

{u, Lixau) < C{E[u] + ||u(l)|P) (7.52) 

□ 

This has proven an estimate for {u,L^Xau)- In the later phase space analysis, it will be necessary to 
estimate mixed derivative norms involving and functions of L and p* . The following corollary is a result 

of theorem l7.6l It controls the expectation value of localised terms. This could have been proven directly 
from the Morawetz estimate with out angular modulation. 

Corollary 7.7. There is a positive constant C such that for all u E S which satisfy the linear wave equation, 

' 1 

{u\ Y—^u')dr < C{E{u) + \\u{l)f). (7.53) 

In particular, if u is a solution of the linear wave equation, il + Hu — 0, and xi '^s a smooth, compactly 
supported function, then there is a C such that, 

{u\xiu')dT <C{E{u) + \\u{l)f), (7.54) 
\^XiufdT<C{E{u) + \\u{l)r). (7.55) 

Proof. Using the first term in the first estimate of theorem 17.411 with m = gives this result. □ 

We also have control of the full two angular derivatives in any compactly supported region, with a weight 
which vanishes quadratically at the photon sphere, = 0. 

Corollary 7.8. For any continuous, compactly supported function x, there is a positive constant C such 
that for all u Cz S which satisfy the linear wave equation, il + Hu = 0, 

{u,plxL^u)dT < C{E[u] + \\u{l)f). (7.56) 

Proof. Again using the first estimate of theorem l7.41l with m — 0, this time using the second term, and noting 
that any continuous, compactly supported function which vanishes quadratically at the photon sphere, in 
particular plx^ is dominated by —gV[. On this compactly supported region, the difference between — I 
and is order one, so it's a local decay term, and also space time integrable. □ 
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8 Phase Space Analysis 



The main conclusion of this section is the phase space induction theorem. It does not control the time 
integral of the angular energy, but instead loses an arbitrarily small power of L. The bounds are still in 
terms of the initial energy and the norm. 

To prove this result, we introduce the operators Tn.m to majorate . Each r„_m is a phase space 

localised version of 7^™ ■ Our localisation will use the radial variable and the radial derivative rescaled using 
the powers of L. We call these rescaled quantities the phase space variables. Because of the localisation, the 

^ + "^ 

commutator of Tn,m and H dominates L 2 only in regions of phase space. This is where our definition 
of propagation observable differs from the standard one. In the standard definition, the commutator of a 
propagation observable must dominate the operator it is said to majorate; whereas, we allow this domination 

to occur in a region of phase space. Under our definition, but not the standard one, r„^m majorates L 2 . 

The definition of majoration always permits "lower order terms". For r„^m, these will be . So that 

the lower order terms truly are lower order, we require < n < m < ^ and n < j. By combining various m 
and n choices, we eventually show that ||i^^^Xa'"lP is integrable in time. 

In subsection 11.11 Structure of the Paper, there is already an outline of this section, which we repeat 
in brief. Subsection 18.11 introduces the phase space variables and r„^„i. Subsection 18.21 introduces the 
commutator expansion theorem, which allows is to expand commutators of functions of the phase space 
variables, and several lemmas for common cases. Subsection 18.31 computes the commutators to show Tn^m 
dominates powers of L with localisation in both phase space variables. Subsection 18.41 combines these 
estimates to prove estimates which are localised in the rescaled radial derivative alone. Finally, subsection 
18.51 combines these in a finite induction to eliminate all phase space localisation. This is sufficient to prove 
pointwise in time decay. 



8.1 Phase Space Variables, Localisation, and Multipliers 

Previously, we introduced a radial variable which was rescaled by the angular momentum. Now, we also 
introduce a rescaled radial derivative. We refer to these as the phase space variables and respectively. 
These will be rescaled separately with different parameters m and n. By varying these parameters, we will 
be able to control powers of L in various regions of phase space. 

Definition 8.1. The phase space variables are 

x.„ =£"p. (8.1) 

CoroUarv 17.81 provides control on the angular energy away from = 0, so that the angular component 
of the energy only needs to be controlled in a compact set around ~ 0. This is the region near r — a. We 
introduce the weight to work in this region; although, the particular choice of Xa will not matter. 

Definition 8.2. The function Xa is defined to he a smooth, compactly supported, radial function, such that 
Xq !i 1; Xa ~ ^ in a neighbourhood of — 0. 

For the phase space variables, we have a notion of lower order, involving our choice of Xa- 

Notation 8.3. The notation 0{L^) denotes functions of time which are integrable, and for which the 
norm in time is bounded by a combination of constants, the energy of u, and its initial L^(9Jt) norm. By 
theorem \6.b\ corollary \7//\ and corollary \7.t\ this includes local decay norms, localised radial derivatives, and 
angular derivatives with vanishing weight at — 0. 

The notation B denotes an arbitrary bounded operator, just as the notation C denotes an arbitrary 
constant. As with arbitrary constants, the value of B may vary from line to line in an argument. The 
notation Bi is used to refer to a bounded operator which is referred to later in the same argument. 

For n < i, the notation 0{\\L^~ Xau\\^ , L^) denotes inner products and norms of u which are either of 
the following form 

||L^Xau|p,0(ii), (8.3) 
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or bounded by norms of this form. Note that for q < 1 + 2n, this includes any term of the form 

WLiBxcuf. (8.4) 
Note that in particular if m < i, then m + n < ^ + n and 

WL^'+^Bxo.uf - 0(||L^Xawf,i'). (8.5) 

To begin with, it is not possible to use sharp cut off functions in x„i and so more regular functions 
with infinite support must be used. The rate of decay for these functions is important, and various sec- 
ondary functions derived from the also appear. The localisations occur in families, with X referring to x„i 
localisation, $ referring to localisation in a wide band around zero, and ^! referring to ^„ localisation in 
a narrow band away from 0. 

In our estimate on the angular energy, we lose powers of L from our choice of localisation. We lose a 

\x\<l- 



factor of L'^ from our choice of decay for ^a.e{x) and a factor of from our choice of width for '^l-s^i 



Definition 8.4. The weight g is defined by 



bx -j^ 



9i^) = I , l_l^2 ^^ (8-6) 



/o (1 + M) 

where b is chosen as in lemma \7~^ The near and far localisation functions for x™ are 

b 

il + b\x\) 



Xiix)^ ,^ , (8.7) 



1 



X,(,)^,^9M. (8.9) 
From this, X|(xm), X|(xm), and X|(xm) are defined by the spectral theorem. 

Definition 8.5. The smooth, near localisation for ^„ and two functions derived from it are defined by 

$a,,(a:) =(l + a;2)-^ (8.10) 
$fc,,(x) ^x^'U^) (8.11) 



*c,.(a:) =y'$a,e(a:)($a,e(a^)+2a:$^_,(x)). (8.12) 
(The definition of ^c.e{x) requires the argument of the square root to be positive, which is shown in lemma 

We use the notation x(^, x) to denote the characteristic function of A which is 1 if x € A and otherwise. 
The sharp, near localisation and sharp, interval localisation for 4„ are 

$|.|<i -X([0,l],x) (8.13) 
*/--<|.|<i =x(r',l],2:). (8.14) 

Smooth, interval localisations will be defined in section \E~^ The functions ^^x\<i '^'^c? ^'/-'S<|2:|<i extended 
as even functions. 

From these, 'l'a.£(^n), ^6.e(^n), $c.c(^„), ^|a;|<i o^*^ ^L-*<|^ |<i Can be defined by the spectral theorem. 

To prove a phase space localised version of the Morawetz estimate, we introduce the operator Tn,m- The 
parameters n and m are restricted to < n < m < i , because this is the range on which an estimate can be 
proven. 

Definition 8.6. For < n < to < i, the phase space induction multiplier is defined to be 

rn,m =Xa<i>a,.(|„)i"~'7L™'i>a,e(|JXa (8.15) 
= ^Xa$a,.(eJ« (5(X™)|„ + ^„5(X„0) <^aAin)Xa (8.16) 
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We now prove some preliminary results, starting with some estimates on the localisations. 
Lemma 8.7. For all e > 0, n e [0, \], and w g § 

\\^c,Mn)v\\ >eH<^aAU^l (8.17) 

W-^ir.'^aAUM <\\L-^^V\\ + II A^ll (8.18) 

Proof. The function '^c,t{x) is computed explicitly. 

^,,,{xf =^a.AUi^aAU + '^X^'aA^)) (8.19) 

= (1 + x^)-^ (1 + ((1 + a;2) - 2^-j^2x'') (8.20) 

^{l + x^y^ {I + €x^) (8.21) 
>e{l + x^)-^ (8.22) 
=e*a,.(0' (8.23) 

Since <i>a,e($„)($a,e(Cn) + ^a;^^ ^(x)) is strictly positive, the square root function, ^cA^) is well defined. 

Since <^cA^Y ^ ^^aA^nY^ the spectral theorem also determines that ||$c,e(€n)^l! ^ ll*J'a,£(^„)?^||. 

Equation H8.18(l is proven first for Schwartz class data. For this data, in a representation where L and 
— both act as multiphcation operators, $a,e(^„)^ < l^rJ^^^'^ ~ I — ^af^-^"^^!^^^"^. This can be used to 

dp. 



start estimating the norm of L" w—^aA^n)'^'^- This gives equation (|8.18|) . 



\\L'-^^n^aAUM <\\L'-^Uta\-'^'A\ (8.24) 

<l|i'-1l„r^ll (8.25) 
<||^i+„.-2e| _ . 9 i^^ll (8.26) 

Continuing in a representation where both L and —i-^ act as multiplication operators, the identity, for 
positive / and g, fg < + g» can be applied with / = 7^i+"e-2« g^j^^j g = \ — i-^l"^. 

\\L'-^^r.<^.AUM <\\l'^^v\\ + W^^vW (8.27) 

□ 

Since Tn,m is a product of operators, there is a Leibniz (or product) rule for computing the commutator 
with another operator. It is possible to further simplify this by expressing a sum of complex conjugates as 
a real part. The main application will be when the other operator is H or one of the subterms Hi. 

Lemma 8.8. For 0<n<m<i, and any self-adjoint operator G which commutes with L 

(u,[G,r„,„» ^{u,Xc.^aAin)[G.lLAL''^"^aAin)Xc.u) (8.28) 

+ {u,xAG,^aAUhL^-L''~'^aAUXc.u) (8.29) 

+ (u,Xa<i>a,.(l„)7L'"i""'[G,<i>a,e(IJ]XaU) (8.30) 
+ (u,Xa<i>a,e(l„)7L'"i""'$a,.(l„)[G,Xa]ti) (8.31) 
+ (ti, [G, Xc.]^aAUlL^-L^''^aAin)Xo.u) (8.32) 

Proof. This is simply an application of the Leibniz rule for commutators. □ 
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8.2 Commutator Expansions 

To apply the Heisenberg like relation, it is necessary to expand commutators involving localisation in 
This is done through a version of the commutator expansion lemma previously used in scattering theory]^. 
We consider the special case of this expansion for commutators involving localisation in the phase space 
variables ^„ and or Xm- These expansions are as finite order power series with an error term which 
involves the Fourier transform of a fcth order derivative of one of the localising functions. 

First, the adjoint is defined. These are iterated commutators and will appear in the commutator expan- 
sion, both as terms in the finite power series expansion and as a term in the remainder. 

Definition 8.9. For two operators H and A, the kth commutator of A with respect to H , Ad\{H) is defined 
recursively by 

Ad\{H)^[H,A] (8.33) 
Ad'x{H) = [Ad'r\H),A] (8.34) 

Initially the commutator A] is defined only as a form on the domain of A intersect the domain of H by 
the formula 

{u, [H, A]u) ={Hu, An) - {Au, Hu) (8.35) 

If Ad'^^{H) extends to a bounded operator, then Ad!\[H) is defined on the domain of A. 

The commutator expansion theorem expands the commutator of an operator, Fi , with a function of a 
self-adjoint operator, F2{A). The expansion is as a power series in the adjoint and has remainder involving 

1 . \k] 

the L norm of Fourier transform of ■ 

Theorem 8.10 (Commutator Expansion Theorem). If n > is an integer, A is a self-adjoint operator, 
Fi is a self-adjoint operator satisfying 

1. for 1 < k < n, Ad\{Fi) extends to a bounded operator, 
and F2 (x) is a smooth function satisfying 

2. |l.F[i^["l]||i<oo, 

then if [Fi^ F2{A)] is defined as a form on the domain of A^^ , 

n—l ^ 

[F,,F2{A)] = J2 ^Ffk^)^4(^i) + Rn (8.36) 
fc=i 

in the form sense with the remainder i?„ satisfying 

\\Rn\\ <C\\:F[Ft%\\Ad^A{Fi)\\ (8.37) 
Consequently, [Fi, F2{A)] defines an operator on the domain of A"^^ . 

Proof. This is proven in [221 ■ '-' 

The commutator expansion theorem is now specialised to the localising functions in the phase space 
variables. In terms of the previous theorem, A = = —i-^L"'~^ and Fi is a function of or of x.,„. 

There are two applications of this lemma. The first application is the expansion of the commutator of 
the Hamiltonian, H, and the phase space induction multiplier, r„^m, as a power series in the adjoint. The 
remainder is a higher order term, in the sense that the commutator minus the power series expansion can be 
multiplied by a positive power of L and still remains a bounded operator. The second application is using 
the expansion with no terms and using the remainder to directly estimate the entire commutator as a higher 
order term. This second application is expanded upon in the lemmas which follow this one. 
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Lemma 8.11 (Commutator order reduction lemma). Let k be a positive integer, Fi{x) be a smooth 
function satisfying 

1. for l<j<k, \\f[^^\\oo < oo, 

and F2 (x) be a smooth function satisfying 

2. forl<j <k - 1, ||fJ^''||oo < 00, 

3. <oo 

then there is a constant, Ck, depending only on k, such that 

k-l 

L'^^'-'^^FM.F.iU] = ^L'=(^-"'Fi^'(|J(iL"-i)^Fp'l(p,) + (8.38) 

i=i 

<C..||.F[Ffl]||i||FW|U (8.39) 



fe-i 

i'=('-™-")[Fi(x,„),F2(^„)] = ^L'=(i-'"-")Fi^l(|„)(zL"+"-i)^Fp'l(p,) + Rk (8.40) 

||i?fc|| <Cfc||^[Ffl]||i||i^fi||eo (8.41) 

Proof. The main idea in this proof is to apply the commutator expansion theorem, theorem 18. 101 
The m = case is a special case of the formula for general to > 0. 

All quantities in equation H8.41|l are composed of powers of L and functions of operators which commute 
with L. Therefore, they preserve the spherical harmonic decomposition, and it is sufficient to prove the 
power series expansion of the commutator on each spherical harmonic. We will use I to denote the action of 
L on each spherical harmonic, rather than the standard spherical harmonic index. 

On a fixed spherical harmonic, take x — = l"^p* and A = = and apply the commutator 

expansion theorem. 



.=1 

\\Rk\\ <Cfe||.F[^^fl]||i||Ad^^jFi)|U 
The adjoint can be rewritten in terms of I, p^, and 

Adl (Fi (x„. ) ) = [Ad^^;; ' {Fi (x„) ) , U 
= [Adl:\Fi{l^p^),~il 



k-l I 

(8.42) 



fc-i ^ 

Fi(x„0,F2(|J] -Ff l(|JAd|jFi(x,„)) + R, 



Tn-i d 

dp.. 



dp* ^" 



Taking a derivative with respect to of a function of l"^p. will introduce an additional factor of U 
j = 1, we have 

Ad^ (Fi(x„0) =^l"-'^F^{^p,) 
^" dp, 

=if'"+"-iF{(r™p,) 

Each adjoint acts as a derivative, so, by induction, 

Ad|^ (Fi (x„o) =iir+"-yFl'^ (rp.) 



8.43) 

8.44) 
8.45) 

8.46) 

. For 

8.47) 
8.48) 

8.49) 
8.50) 
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Applying this, and multiplying by l'^^^ " gives that on each spherical harmonic, 

f'^™+"-'H^^i(x™),F2(ej] = ^ -Ffl (In) (x™) + (8.51) 

\\Rk\\ <C,||^[i^fl]||i||i^[^'l|U (8.52) 

Since the constant in the commutator expansion theorem, theorem 18. 101 is independent of the choice of x, 
A, Fi, and F2, the constant Ck is independent of I. The result on each spherical harmonic can be extended 
across all harmonics with out any change. In particular, the operator R^. is uniformly bounded across all 
spherical harmonics. □ 

The commutator order reduction lemma, lemma IS. Ill is used to estimate commutators of localisation in 
the phase space variables. Usually it is sufficient to know that the commutator is a bounded operator times 
negative powers of L, and that is what the following lemmas assert for certain classes of localising operators. 

The following lemma applies in the common situation when the localisation in ^„ is a Schwartz class 
function. 

Lemma 8.12. There is a constant C, such that if < n^, < m < i, Fi G C^(M), F{ G , and F2 is of 
Schwartz class, then there is a bounded operator, B, satisfying the following 

[i^i(x„),i^2(IJ]=L"+"-^Si (8.53) 
Pill <q|.^[i^:;]||i||i^i'||oo. (8.54) 

Proof. Since the Schwartz class is preserved by differentiation and the Fourier transform, .^[^"2] Schwartz 
class and, hence, in L^. The expansion to order from the commutator order reduction lemma, lemma lS.llI 
is 

Li-'"-"[i^i(x„),F2(|J]=i?i (8.55) 
||i?i|| <C||^[F^]||Li||i^:;i|L~ (8.56) 

Since L is bounded below by 1, and < n, m < i, the operator is bounded. The operators 27™+"^^ 

and i?i are bounded, so they can be composed. Applying 2^™+"^^ to H8.55|l proves the desired result. □ 

The previous lemma can not be applied with the smooth, near localisations for which are not in 
Schwartz class. To prove a similar result for the smooth, near localisations, it is sufficient to show that each 
satisfies the condition J-[F2] G L^. The following lemma also proves a commutator order reduction result for 
two unbounded functions of F2 — S. and F2 = ^$o,£(0' using a first order expansion in lemma IS. Ill 

Lemma 8.13. There is a constant C , such that if < n < ^, < m < ^, Fi has its first in , e > 0, and 
F2{x) is one of the following functions: x,<^a,e{x),<^}j^^{x),^c,t{x), then there is a hounded operator, Bp^, 
satisfying the following 

[i^i(x„0,F2(^„)] =U^+'"^-'Bp, (8.57) 
IIBf.II <C\\F[\\^. (8.58) 

// the previous conditions apply, Fi also has its second derivative in L°° , and F2{x) = x^a,eix), then 
there is a bounded operator, Bp^, satisfying 

[Fi(x„),F2(|J] ^U'+^^'-^Bp, (8.59) 
\\Bp^\ <C{\\F[\\^ + m\^). (8.60) 

Proof. The fundamental problem here is that the permitted F2 functions are not in Schwartz class. There 
are finitely many functions F2 considered in this theorem, so if a C can be found for each of them, then the 
largest one can be applied to make the result hold uniformly. 
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For the function F2{x) — x, i^2(^n) = = ^*gf~^"~^- operators involved commute with 

the spherical harmonic decomposition, and on each spherical harmonic, the commutator can be explicitly 
computed 

[i^i(x„), -^^L"-^] ^[FiiL'^p,), -*7r-^"~'] (8-61) 
op* op* 

=iL"+"-ii^^(L"V,) = iL"+™-iFi'(x„.) (8.62) 

The commutator is a product of a power of L and a multiplication operator. The L°° norm of the function 
is the operator norm of the associated multiplication operator. 

To apply the commutator order reduction lemma, lemma 11^.111 to the remaining functions, we prove: If 
h e C2(M) n L^(R) and h" e L^{R), then \\T[h]\\i < oo. 

Since h ^ C^, h" is well defined, and since h,h" E L^(M), both J-[h] and are well defined and in 

L^(R). Let 

f{o^\{i+emhmi (s.es) 

/big(0^/(Ox([l,oo),/(e)), (8.64) 

/sn.an(e)=/(e)x([0,l),/(e)), (8.65) 

and observe that 

/(0-/big(e)+ /small (0 (8.66) 

\mm <\^\ + < i/big(oi + (8.67) 



From properties of the Fourier transform it is known that 

\nhm\ + mh"m\ = imm + \emm - 1/(01, (8.68) 

and hence that / is in L^{R). Since /big < /^jg < /^, /big is in _L-'^(R), and since /small < I7 '^""^'gi^'* is also 
in L^. The sum of these bounds so that J-[h] is in L^(M). 

For the three functions, $a,e(a;), ^b,eix), and ^c,e{x), each of these is a power of a polynomial and decays 
like x~2— . Therefore, and (F2)" ^''^^ each in L^, since they decay like x~^ and x^^ respectively, and 
^^[Fj] is in L^. By the same argument as in lemma IH. 121 

[Fi(x,^),F2(|J]=L«+™-iB (8.69) 
\\B\\ <C{\\Fi\U. (8.70) 

The last, and most difficult piece, is for F2{x) — x$a,e(a;). In the previous arguments, the commutator 
was directly estimated. In this case it is expanded to first order, and both the first order term and the 
remainder are estimated. The derivative of x^a,e{x) decays too slowly to be in and may not have an 
Fourier transform. 

As usual, all the operators involved commute with the spherical harmonic decomposition. On each 
spherical harmonic, by lemma lS.llI 

[Fi(x,„),F2(CJ] =F^(CJ(zL"+™"i)(F{(x,„)) +i2"+2™-2^^^ (8 

The first term is a product of operators which commute with powers of L and each can be estimated in 
norm. 

mui-^L-+"^~')i-Fii^^m =L-+^-'mu\mM\\ (8.72) 

<i«+™-i||F^||^||Fi'|U (8.73) 

Since ^2(2;) = x^a.e{x), F2 is bounded. 
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The second term is estimated by lemma IS. Ill 

\\R2\\^C\\:F[F^']UFi'\U (8.74) 

Since F2 — x^a,e{x), F2 is a smooth function, which is in and has its second derivative, {F^)" , also in 
L^, we conclude that .^[Fj'] is in L^. 

Since x2n+2m-2 j^n+m-i ^ ^^iQ estimate on each of the terms in the expression for the commutator can 



be combined to give 

||[Fi(x„0,i^2(IJ]|| <L"+™-^||F^||oo||Fi'IU (8.75) 

<CL^+"'-W\F[\\^ + \\F^\U) (8.76) 

On each spherical harmonic there is a bounded operator, B, such that 

[Fi(x„),F2(^J]=L"+"-ii? (8.77) 

||S|| <C{\\F[\\^ + \\F^U). (8.78) 

Therefore, B can be extended as an operator on LP' satisfying the same condition. □ 



8.3 Phase Space Estimates 

In this subsection, we compute the commutators between Tn^m and the components of the Hamiltonian, 
Hi. This is analogous to subsection 16 . 21 with phase space localisation and additional factors of L. The H2 
commutator is simply lower order. The Hi + H^ commutator dominates the radial derivative localised to 
small and powers of L localised to large x^. 

The commutator with H2 is computed first. Because H2 is a function of and contains no derivatives, 
all the commutators can be shown to be lower order by lemma IS. 131 

Lemma 8.14. For u E §> is a solution to the wave equation, il + Hu = 0, and < rt < to < 

(u, [H2, ) = 0{L^) (8.79) 

Proof. Lemma 18.81 the Leibniz formula for commutators, is used to calculate [-ff2,r„^m] and then the com- 
mutators are shown to be lower order by lemma 15. 131 Since Xa and V commute, [i?2, Xa] = and only three 
terms are left in the expansion of the [-ff2,r„^m]. 

(U, [i/2, r„.™]u) = (U, Xa L^-\[V, ^aMn)]lL'^^^aAU (8-80) 

+ ^aAU[V,lLA'^aAU (8.81) 
+ ^aAUlL^ ^a.,Mn)\)Xc.u) (8.82) 

The commutator [V, 7i,m] only involves the angularly modulated multiplier, and was shown to be a bounded 
operator in lemma 1731 Applying additional bounded operators and powers of L yields a bounded operator 
multiplied by the same power of L. 

L''-''^a.Mn)[V,lLrr^]-^a..Mn) = L"~'B (8.83) 

Since n < § , by the angular modulation theorem, theorem 17.61 the expectation value of this operator is 
0(Li). 

{u,Xc. U'-''^aAU[V,lL-]'^aAUXaU) ={u,V'-'Bu) (8.84) 

=0(Li). (8.85) 

The remaining terms are complex conjugates of each other. Only one is considered here, and the other 
can be dealt with in the same way. The commutator is shown to be a bounded operator times a negative 
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power of L using lemma IS. 131 71,™ can be expanded and the bounded functions involved in this expansion 
can be absorbed into the bounded operator from the commutator expansion. 

L-'-^aAUlLr^ [V. -^aAU] ^L^-'^aAtn){^jL^ " ^5^™)^""'^? (8-86) 

=L2"-i-^ci>^^^(|J_^Bi + L2"-i-^L™ixx(x„)i?2 (8.87) 
op* 

This operator can be substituted into the relevant expectation value. This gives a sum of two terms, both 
of which can be estimated by the Cauchy-Schwartz inequality. 

{u, Xc^L^-'^aAUlL^ [V, ^aAU]Xc.u) (8.89) 

H^^^aAUXcU, L'^-'-'BiXcu) + ixcu, L^+^^-i-^Bxaw) (8.90) 
op* 

- ~ (^(Wq^XcuUL^^BiXcuW + ||L^"xau||2) (8.91) 

- ' ^^"9^^"""' + ll^'^Xc^lP + WL-^'^XcuW') (8.92) 

Since 2n — 1 — e < | and (3 — e)n < f , by the angular modulation theorem, theorem 17.61 the two terms 
involving powers of L are time integrable, ie 0{L^). By corollarv l7.7l Hg^XaulP = O(L^). Therefore the 
expectation value being considered is time integrable 

{u,Xc.L''-''^aAUlLAV,<^aAU]Xc.u) =0{L^) (8.93) 

□ 

Because = (— A52)Vl = (L^ — 1)Vl contains two angular derivatives, which are exactly what is under 
consideration, this commutator requires the most care. The difference between — 1 and is a local decay 
term. The commutator is first expanded using the product rule in lemma 15^ The commutator with ^^m. 
is computed from the angular modulation argument in lemma 17^ The commutator with the localisation 
'^aA^n) can be expanded to second order using the commutator order reduction lemma, lemma IS. Ill The 
first order term is of the same order as the commutator with 7l">. It involves a derivative of <^aAin) and 
is multiplied by the term from 7^™. There are two such terms, and when rearranged and combined 

with the term from the commutator with 7l™, they give $c,e(^„)^ = *I'a,e(€n)(^a,£(€n) + 2*I>a,e(€n))- This 
term must be positive to give the desired estimate. This dictates the choice of the decay function ^aA^) 
and forces e > 0. The second order term in the expansion naively appears to be of a high order since 7/3 
introduced two factors of L. However, rearrangements yield a cancellation of the highest order part of this 
term and the remainder is a lower order error term. For both the first and second commutators, there are 
many rearrangements of phase space localising functions. These are shown to be lower order in the standard 
way using lemma IS. 131 The remainder term is shown to be lower order directly from the commutator order 
reduction lemma, lemma IS. Ill 

Proposition 8.15. For u £ S is a solution to the wave equation, il + Hu = and < rt < m < 

{U, [Hi + iJ3,r„,„» >c($a,e(€„)XaW, (-^X|(x„0^ " ei'sL" I^L )*a,. (4„)XaU) + 0(||L^XaW||', L^) 

(8.94) 

Proof. The commutator is given by the Leibniz rule as the sum of five terms. 

+i?3,r„.™] ^H, + H3,Xo.]'^aAUL''-'lL^^aAUXc. + Xa^aAUL'^^'lL'^'^aAUlHl + H3,Xo.] 

(8.95) 

+ xAHl + H3, 'faAU]L"~'lL^^aAUXc + Xa'^'aAUL'^-'lLAHl + i?3, $a,e(l„)]Xo 

(8.96) 

+ Xa$a,.(l„)i"-Mi^l +i?3,7L">]$a,.(€„)Xa (8.97) 
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Since the first two terms are adjoints of each other, it is sufficient to estimate oniy one of them. The 
remaining three terms are more compiicated and must be estimated together. 

Since Xa is a function, it commutes with Vl, and we only need to determine the contribution from —^3-, 
which can be expUcitly calculated. 

(w,Xa<fa.e(l„)7L'"i""'<i>a,.(€„)h^,Xa]«) (8.98) 

= (u, x„$,_,(|J^^„L»-^cl>,,,(| J(-2 A;^'^ _ (8.99) 

op* 

To analyse this term, is expanded into two terms so that the most significant term can be dealt with 
first. The slightly less common expansion 7^™ = -^gL^ — ^g'^m is used. There are a total of four terms to 
consider, two from expanding [Hi, Xa] times two from expanding 7^™. 

The term that appear to be highest order is the one involving -^gL"^ from "fi^m and ~2^^x'a from 

{—-^,Xq\- It will be referred to as /q. It can be simplified by moving one of the derivative operators 
through <^a,ii{i,n) (since they commute), and then moving it and a factor of Xa to the other side of the inner 
product. 

/o ^(^i,x.$a,.(IJ^.9L".i^"-'$a.e(l„)(-2^x:.)") (8-100) 

= - i^XcU, ci>,,,(^Jgi„L"-^$„,,(^ J(-2 A^^)^) (8.101) 
ap* op. 

Applying the Cauchy-Schwartz estimate gives the sum of two terms. One of these is llg^X'^ll^ which is 
0{L^) (integrable in time) by corollarv l7.7l The other can be rewritten in terms of 

h>- W-^Xcuf ~ ||$,,4|Jgi„L"-^$,,,(|J A;^'^„||2 (8.102) 

> - ||il-^$a.e(IJ.9L™l„$a..(IJX>ll' + 0{L^) (8.103) 

The remaining norm is going to be manipulated with the goal of applying lemma I^TI to control ^„'I'a,e(^n)^ 
by powers of L. This manipulation will introduce additional error terms from commuting the phase space 
localisations. To commute *l'a,e(^„) and while avoiding any error terms involving factors of ^„ outside a 
commutator, it is necessary to commute ^„ through 51,™ and then commute ^„^a,e(^„) back through gh^. 

Io>- \\L'-%,<^aAU9L^'i'aAUx>f (8.104) 

- ||Li-^$<,^,(|J[.9L..,e„]$a,.(€Jx>ll' + 0(i') (8.105) 

> - \\L'-'gL'^^n'faAU\>\\^ - \\L'-ntn'faAU,9Lr^]'^aAtn)X>f (8-106) 

- \\L'-'<faAU[9L'^,U^aAUx'c.u\\' + 0iL') (8.107) 

By lemma 1^.131 the commutators are of the form L"'+"~^i3. Since <I>a,e(^„) is also a bounded operator, all 
the error terms from commuting are lower order terms. In the remaining term, the factor of can be 
dropped for a constant and then lemma IHTI can be applied to control ^„$a,e(^„)^- 

Io>- C||Li-^5i™|„$,^,(^ J2x>||2 - C||Li-^i"+™-ii?x>||' + Oil') (8.108) 

> - CiWL'^^x'c.ur + II^X>ir) - OiWL'^XcuW^L') (8.109) 

The angular energy away from r = a is time integrable by corollarv l7.8l Since Xa is compactly supported 
and identically 1 in a neighbourhood of r = a, |XqP is dominated by a multiple of {p*g{p*)) and HLXa""!! is 
integrable in time. Interpolating between the local decay result and the angular energy decay away from the 
photon sphere result, \\L r^^^ xL^lP is time integrable, and hence in 0{L^). By corollarv l7.7l Hg^xL^lP 
is time integrable and hence in 0{L^). All the terms in /q are now controlled. 

/O ^(z,,x.<i>.,e(CJ^5L™L"-^<i>a,e(€J(-2^x'J") > 0( 1 1 Xa^ll ' , ) ■ (8.110) 

ap* ap* 
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The term that appears to be the next highest order is also 0{\\L^^ XauW^ , L^) not merely 0{L^). This 
is the term involving from 7/,™ and —2-^x'a from {—-§^,Xa\- estimated by moving most of 

the localisation to the other side of the inner product, applying the Cauchy-Schwartz inequality, replacing 



bounded operators by a constant, and then finding that the remaining terms are each 0{\\L ^2 Xq^^IPi L^) 



1 L"+"-^ j)_ , 

2(1 + L™|p,|)2^ 



{u, X«'i>a,.(€„)(-777^^;;^^-^)<i>a,e(^J — X» (8.111) 



= - i(<i>a,.(|Ji™+"-'Xi(x™)$,,,(CJxa", (8.112) 



>-C||L'"+"~^Xau||||T^X>|| (8.113) 

> C{\\L"^+"-^Xaur + II T^X>|P) (8.114) 

op* 

>-Oi\\L^Xo.uf,L') (8.115) 

The terms involving Xa from [—-^jXa] £^re the lowest order terms and are quickly shown to be lower 
order. The same method as was used for the term involving and —2-^x'a can be applied. 

dp* op* 

>-C\\^Xc.u\\\\L--\y\\ (8.117) 
op* 



>-0{L^) (8.118) 
=0{\\L-^Xc.u\\\L^) (8.119) 



jjm-\-n — e 

{u,Xo.^a.AU-—^^L^^aAUX» (8.120) 

H^a.AUXcU, L™+"-^Xx(x„0*i>a,e(^ JX» (8.121) 
> - C\\L'^^aAUXo.u\\\\L'^^aAUXo.u\\ (8.122) 

>-C||L''^Xc."||' (8.123) 
>~0{\\L-^Xo.u\\\L') (8.124) 

We now turn to the three remaining terms in (|8.95|l . For the commutators with $a,e(Cn); since —-§p! 
commutes with it is sufficient to compute the commutator with Vl- No distinction is made between 
[i73,r„^m] = \VAL'^ — 1),!"™,™] and \VlL? ,Tn,m] since the difference is lower order by an identical argument 
to that in lemma IS. 141 These commutators are expanded to third order. 

[VL^'^aAU] = <Mn)[VL,U+KAUWL,U,U+R:^- (8.125) 



In each of the third and fourth terms of (|8.95|) . we make this substitution and expand 7^™ as a sum of two 

terms. This gives twelve terms, which are combined with the fifth term in (|8.95|l . for a total of thirteen. 
These are grouped as 

XAHz^'^aAUWlL^'^aA^UXo. (8.126) 

+Xa<i>a,e(^„)i"-^7L'" [^3, ^a,. (I „ )]Xa (8.127) 

+Xo.^aAUL''-'[Hl+H-i,^L^]<^aAUXo. =Xa<i>a,e (IJi""^ h ^ , 7L'"] I'a,. (OXa (8.128) 

+ Xa(/l +h+h+h+h+ h)Xo.- (8.129) 
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where 



h 



dp* 



d 



op* 

d 

= L^+"-'{^'UU[[VL,U,U^jL'^^^aAU 

d 

+ ^aAU9L^-g^^KA^nWL,U,U) 

L'+^-'\{ KAU[[VL,U,U9'L^^^aAU 



+ ^aAU9'L^KMn)[[VL.UXn]) 

op* op* 



(8.130) 

(8.131) 

(8.132) 

(8.133) 

(8.134) 

(8.135) 

(8.136) 
(8.137) 
(8.138) 

(8.139) 



The leading order term^, /i, is composed of two types of terms. One type comes from the commutator 
[L^VliTl™], and the other comes from combining the first order term in the expansion of [H^,(^aAin)\ 
with the derivative terms in 7^™. To get the desired estimate, it is necessary to rearrange the phase space 
locahsing functions, which introduces additional lower order error terms. In the calculations, leading order 
terms are written first here, followed by those which will be shown to be lower order. 

The term Ii is rearranged so that all the leading order subterms contain the p* localisation factor, 
(— V^(7im), which first appears in the commutator [Vl,7l™], and is the localisation in the statement of the 
lemma. The commutators from rearranging localisations are lower order by lemma 15. 131 



2+n-€ 



{^aAU{-Vi){-^L'^-')4-9L^^aAU 



dp* 



d 



+ '^aAU9Lr^-^^'^aAU{-vL){-^L--A) 

= L^^''-%<^aAUU'VigLr^)<^,AU 

+ ^aAU{-VLgL^)^aAin) 

+ -^aAU(~VigLr<^)^^'^aAU) 

+ L^+--^{<^aAU[-Vi,U9L^^^^aAU + 
<^aAU9LAin-^aAU,-VL]) 

= L^+-'{i^^,AU{-vLgL^)^aAin) 

+ ^aAUi-VLgL^^)<^aAU 

+ ^aAUi-VLgLr<^)^n-^aAU) 



8.140) 
8.141) 

8.142) 

8.143) 
8.144) 
8.145) 
8.146) 
8.147) 
8.148) 
8.149) 
8.150) 
8.151) 



Following this, the ^„ localising functions are gathered to the right of the space localisation. Once they 
are gathered this way, their sum is the operator $c,e(^„)^- One factor of ^cA^n) <^an be moved back through 

®From here to the end of the estimate on /g, the proof is cut and pasted from the proof of prop 8.15 in the previous notes. 
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the localisation, so that $c,£(^„) appears on both sides of the p* localisation. Rearranging in this way 
introduces new commutators, which are not immediately estimated since they involve (— V^iffL™)- This is 
a localising function in and so that the previous lemmas which estimate commutators can not be 

applied. 



h = L'+"-'{ - VlgL'^)<^>aAU{'^aAtn) + 2|„<i>a,. „ )') (8.152) 

+i'+"-'(K„$a,e(l„)', -V[gL^]^aAU (8.153) 

+ [$,,,(!„), -F/5L."]$a,e(IJ (8.154) 

+ [^aAU, -V[gLAi„^aAW) (8.155) 

j_j^l+2n-e^ (8.156) 

= L^+'^-'i^cAUi-vLdL'A'^cAU) (8.157) 

i'+""^([(-^/5L-), ^cAU]'^cAU (8.158) 

+ [in'^aAW, -VlgL'-]<^aAtn) (8.159) 

+ [<^>aAU,'V[gLA^aAU (8.160) 

+ ['i>a,.(l„), -l45L"]l„$a,.(^„)') (8.161) 

_^^l+2n-e^ (8.162) 



Of the terms in the previous expression, the first is a positive term we wish to keep, and the expectation value 
of the last is 0{\\L^~ XauW^ , L^)- The other terms all involve a commutator of the form [-F(^„), (— V^ffL'")] 
multiplied by T^^+n-e bounded operators. 

By rewriting the commutator [-F(^„), (— V/^<7l™)], the four remaining commutator terms can also be 
shown to be 0{\\L^^ XauW^ , L^). It is necessary to separate —V[gL"^ into a function of p, and of p^U^ to 
estimate the commutator. Naively applying the Leibniz rule for commutators, gives a contribution from the 
commutator with which is l^+^^+'^^^-^b, which is not 0(||L~2^Xa^^ll^) iii expectation value. The 
decomposition 

-VLgL^^ = {-Vip-A){p.L^gL^)L-^ (8.163) 

ensures that the commutators are all 0{\\L~^ Xau\\^ tL^) in expectation. This decomposition is possible 
since —V[ vanishes linearly at p* = 0. Taylor's theorem guarantees ^V[p^^ is C°°. 

l2+"-[F(|„), i-VlgL^)] =l2-'"+"-[F(|J, i-V[pA)L"'p4 (8.164) 

The Leibniz rule can now be applied and the resulting commutators shown to be lower order operators by 
lemma 15. 131 Note that the class of F considered here consists of ^c.e{x), ^b,e{x), and $a,e(2^), all of which 
are covered by lemma 15. 131 and that (V^p^^) and L™ p^.gLm = X|(x,„)^ both have bounded first and second 
derivative as required by lemma IS. 131 

L^+'-'WiU, (-^/5L"0] ^L^-^+^^-^{-Vlp-A)[F{U^X^{^raf] (8.165) 

+ L^-^+n-.[F{U.{-VipA)]X,{^..f (8.166) 
=L2+"-^L"+™-1S + l2-'"+"-'^L™+"- V*-B (8.167) 

^^l+2n-e^^^^^l+2„-e^ (8.168) 

Since Xa has compact support, and applying the self-adjointness of p*, we are left with an error term like 

\{p,Xo.u,L-^^Bxo.u)\ <C\\L-^^Xo.u\\ (8.169) 
^0{\\L-^Xc.u\\\L^) (8.170) 

This gives the following estimate of the leading order term, Ii . 

{u,Xc.IiXo.u) = \\L'^{~V[pA)^X^{^„^)<^,AUXc.u\\' + 0{\\L'^Xc.u\\^L') (8.171) 
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Using the estimate in equations (|8.164|l - l|8.171|) again, we can rearrange the localisation as 

{u,Xo.hXo.u) = ||L^<i>,,,(^J(-y/p;i)^Xt(x™)xa"ll' + 0(||L^Xauf (8.172) 
The localisation ^cei^n) can be replaced by e^/^$a,£(^„), so the estimate becomes, 

(7/,Xa/iXc,") =e||i^$a,.(C„)(-T^/p;')^XT(x„Oxa"ll' + 0(l|i'^Xow||',i')- (8.173) 
To conclude the estimate on /i, we commute the localisation again, 

(^,X«^iXaw) =e||^^(-V^/p:')^^T(x™)'i>a,.(l„)Xaii|P + 0(||i^Xowf,i')- (8-174) 

It now remains to show that I2, are 0{\\L^2— ^^uW^ , L^). Discussion of these continues from 

highest order to lowest. The first of these is I2 which contains terms involving and ei^n) terms. 
Note that because two different expansions of 71,™ were used in equation H8.129|l . the signs on the two factors 



of g'j^m are different. To simplify I2, the factors of are expanded. 

I2 = L'+''-'^{~KAU[VL,U9'L-'^aAU (8.175) 

+ ^aAU9l^KAU[VL,U) (8.176) 

= L'+"-'l{ KJUHL--'){-VL)L"'X^{^„,)'^aAU (8.177) 

+ <PaAUXli^rn)KAUHL--')hV[)) (8.178) 

- _ KAUvLx^{^,A<i>aAU (8.i79) 

+ ^aAUXl{^m)KAUVL) (8.180) 



This leaves an expression for I2 as the sum of two terms. The goal now is to rearrange these terms so that 
the highest order parts from each cancel. As usual, the rearrangement will introduce commutators. Naively 
attempting to commute ^a,e{^n) o'' ^'a ei^n) with g^m will leave a commutator of the form which 
is not 0(|jL^'i~Xa'i*|Pi ^^) in expectation value. Once again, an additional factor of p*p~^ is introduced 
and used to absorb a factor of L™ into a function of p^.L"^. 

h = '-L'+'-'-{ - KA^Ap.L"^X^{K„.){V[p:')'^aAU (8.181) 
+ $„„(|Ji™X4(x„,X,,(|„)p,(F/p;i)) (8.182) 

Functions of are now moved to obtain the same localising functions, {V[p^^) and p*L™X|(xm) = 
Xm^j,(xm), in both terms. This rearrangement introduces additional commutators which can be estimated 
by lemma 1^.131 

I2 = - KAUi^mX^i^rn))^aAU(VLP*') (8.183) 

+ ^aAU{^mXiiic„,))KAUivLpA)) (8.184) 

+ _ <i>;_^(|J(x„Xi(x„0)[(yZp;l), <i>aAU] (8.185) 

+ ^aAUXl{^m)L"'[^aAU,P*]iVLpA)) (8.186) 

^ _ $;_^(|J(x„X,(x„0)<i>..e(^ J(V^>;^) (8.187) 

+ ^aAUi^mXl{^„,))<P',AUiVLpA)) (8.188) 

+ L'"+3"-^B2 (8.189) 

The ^„ localisation is now grouped to the left of all and p*L™ localisation. The leading order terms are 
both ^aA^A^'a ei^n)i^mXi{xm j){Vlp~^), but with oppositc sigus, SO they cancel. The remaining terms 
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are commutators which are found to be L^+S" <^B in the standard way. Although Xi{x) is not C^, we will 
only need to use xXi{x) which has its second derivative in L°°. 

h = ^Li+2"-(-<,,(|J[(x™X,(x„)),$„,,(^J](F/p;i) (8.190) 

+ *a,e(IJ[(x„X^(x„)),<,,(4„)](y/p;l)) (8.191) 

+L^''''Bi + L™+3"-^B2 (8.192) 
= i'»+3"-e^3 _,_ L^r^-^Bi (8.193) 

Because the powers of L involved are all less than 1 + 2?!, the desired estimate on I2 holds. 

iXaU, hXau) = OiWL^Xcuf, L') (8.194) 

Next It, is considered. This involves a second commutator [[Vl, ^„], which can be explicitly computed 
as (F/')(-i'"^'). 

/3 = L^+^-^{<i>l,{U[[VL,U.U^jL^^'i>aAU (8.195) 

+ $a,. iUdL'^^^K. iUiiVL , U , U) (8.196) 

= L'+-'^{^:AUiVn{-L'"-')^jL^^^aAU (8.197) 

+ ^aAU9Lr^^KAUmi-L''^-')) (8.198) 
op* 

To further simplify this, it is rewritten in terms of ^„ instead of 

h - ^L'+^^-\KMn){V^:)^.n9L^^'^aAU (8.199) 

+ <^aAin)9L^inKMn){Vi:)) (8.200) 

This is rearranged to group ^„ and ^'Ai^A together. Since x^'AiA is smooth and decays like x^~ as 
X oo, x$" J (a;) is a bounded function and ^n*!*" e(^„) is a bounded operator. The rearrangement introduces 
commutators which are estimated in the standard way by lemma IS. 131 

h = L'^'"-%A^C.KAU)Vi:9L^^aAU (8.201) 

+ ^aAU9L'^A^^nKAU)^L) (8-202) 

+A+'-^i^AiU[VA -^U9L^^aAU) (8.203) 
= 7,i+2n-^5^ + L3"-^B2 (8.204) 

Since the powers of i involved are bounded by 1 + 2n, the expectation value of I3 with respect to XqU is 
0(||L^Xa^i||^L^). 

To estimate I4, the commutator [[Vl,^„],^„] and the derivative are explicitly expanded to show I4 
is a product of bounded operators and powers of L. 

L2+n-.l( „ ci>A{U[[VL,U,U9'L'^'^aAU (8.205) 



+ ^aA^n )9A^ i^n Wl, U , U ) (8.206) 

= KAUiVL){-L'-')L"^XA^m)^aAU (8.207) 

+ <i>a,,(CJL™X;(x,„X,,(^ J(F/)(-i2"-2)) (8 208) 

^ ^m.+3n-eg (8.209) 

Since m + 3n — e<l + 2n, the expectation value of I4 with respect to XaU is 0(||L^2— 2^^). 
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Finally the terms Is and /g involving R3 = L^"~^B are considered. These are each the sum of two inner 
products, but by the symmetry and anti symmetry properties of the operators involved, it is sufficient to 
consider just one of the inner products for each. The expectation value of is considered, and the operators 
Xa, ^a,ei^n) are movcd to the left side of the inner product. 

{U, Xc'^aAUL^^^'^'j—gL^^RiXcU) ={-J—^a.Mn)Xo.U, L^+^-' QLr^ RsX^u) (8.210) 

op, op* 

The Cauchy-Schwartz inequality is applied and the resulting norms are found to be of a bounded operator 
acting on g^-Xa" ^-^d on L^"~^~'xctU. By corollarv l7.7l || g^-XaulP is time integrable and hence O(L^) and 

OdlL^xlf , Li). Since 4n - 1 - e < i^, HL^-i-^^^tf is 0(||L^Xa"ll', i')- 

(t.,Xa<&a,e(IJi'+"-^7^3L™i?3Xa^i) >- (||<&a,e(€j/-XaZi||||i'+"-^i'"-'SXa"ll) (8.211) 

op* ap, 

> - maAU^Xaur + WL^+^-'L^^-^Bxauf) (8.212) 
op* 

- " + WL^^'-'-'Xcur) (8.213) 

>-0{\\L^Xau\\^L^) (8.214) 

The term /g is estimated in expectation value by expanding g^™ and finding the commutator is a bounded 
operator times 2^"'"+'"^^^'^. Since 4n + m — 1 — e<l + 2n, this expectation value in 0(||L~'2~Xa'"IPi L^)- 

(w,Xai^'+"+""^5L™-R3Xa«) -(u,Xai'+""^^i(x™)i?3Xa«> (8.215) 

= Xai'+'"+"-^i'"-'i3XaW> (8.216) 
= (u,Xai'"+™"^-^SXa«> (8.217) 
=0{L^) (8.218) 

This term is 0{L^) because 4n + TO— 1 — e< |. This shows the contributions from I2, ■ ■ ■ ,Ie are lower order. 
Combining the results so far, we have that 

{u,[Hi+Hs,T^,^]u) >{u,Xa<^aAUL''-'h-^,lL^]^aAUXo.u) (8.219) 

+ {u, Xo.^aAUL''-'ei-V[gL^A^aAUXo.u) + 0(||L^Xa^ll', i') (8.220) 

>{^{L-^, Xo.)u, U'-'[Hi + eH3, 7L'"]$a,.(Xa)w) + Oi\\L^Xo.uf, L^). (8.221) 

From the angular modulation commutator estimate, lemma 17^ we have the desired result, 

i_d_ 
dp 

+ ce{^aAUXc.u,L^+''-'gL^Vi)^aAUXc. + 0{\\L-^Xo.u\\^L^). (8.223) 

Note that lemma 17^ does not require the full —gL^V^ term, since there's a —(1 — e)gL^Vl term appearing 
on the right. Thus, even though we only have —egL^V[ in equation (j8.221|l . this is sufficient to get positivity 
and control over — (e/2)gL™V^. □ 

8.4 Derivative Bounds 

In this subsection, we remove the localisation on the x„i phase space variable. The resulting estimates are 
localised in ^„ only. Our estimates will be in the regions corresponding to *&|»|<i and ^'l-i5<|,|<i. The ffi'st 

lemma shows that r„ i majorates L 2 ^ and the second, that r„^„ majorates L 2 . 

The first lemma is a quick application of lemma l^.151 with lemma IF^ . 1 41 used to show the H2 contributions 
are lower order. 



U,[Hi+H3,Vn,m]u) >c||L"^Xi(x„0*T^$a,e(|„)XaU|P (8.222) 
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Lemma 8.16. For u £S is a solution to the wave equation, ii + Hu = 0, e > 0, n £ [0, and x compactly 
supported, there is a constant C such that, 

{u,[H,T^^^^\u) >C\\L^^(^\^^^\<^xuf-0{\\L-^Xc.u\\\L^) (8.224) 

Proof. By the results of the angular modulation argument, we have that for any smooth, compactly supported 
function, x, 

- Ax^(x„0 A _ L^g^^Vi >L3™x (8.225) 
op* op* 

Taking x to be smooth, positive, compactly supported, and identically one on the support of Xa a-nd m = i, 
we have, from the result of lemmas [8.151 and 18 . 141 and lemma 17^ 

{u, [H,T^^.]u) >{^,^,[i^^)xo.u,L^-\-—X^{^^)—-eL^gL^.Vi)^aAUXc^^^ 

(8.226) 

>e{^aAUXc.u, Li+''~'x^a,e{UXo.u) + 0{\\L-^ x^uf , Li) (8.227) 

>C\\L^^^^^aAUXc.u\\^ (8.228) 

+ q|i^^^L"-iSxaM||2 + 0(||L^Xa"||',^'). (8.229) 

Since the second term, from [x, $a,e(£„)], is also a lower order term, the desired result holds. □ 

Our next goal is a result with localisation ^ l-^<\$, |<i- We begin by introducing new '5 type localisations, 
which are used as smooth approximations to l-^<\»\<i- 

Definition 8.17. The Junction 'Si : [0,cxd) — + [0, 1] is defined to he a smooth C°° function which has support 
on [^,oo) and is identically one on [l,oo). The function ^2 : [0,oo) —t [0,1] is defined to he a C°° function 
which has support on [0, 2] and is identically 1 on [0, 1]. These are extended as even functions. The extensions 
are also Schwartz class, since the original functions are constant in a neighhourhood of zero. The operator 
^'(L^'',^„) is defined by 

^{L-\U = *i(i'^J*2(l„) (8.230) 

The following lemma allows us to replace one $ type localisation with 'J type localisation and to move 
this replacement through localisation. 

Lemma 8.18. Given 6 and e positive, there is a constant Ci such that for all v (li S and for F2 either ^a,e 

or <i>c,e, 

||*L-<|Cj<iHI <\\^iL-',Uv\\ (8.231) 
Cl||L-*(*(i-^.)F2-l)(^„)«|| <U„v\\ (8.232) 

In addition, for n G [0, and m G [0, there is a constant C2, such that for any differentiable function 
Fi and for F2 either $a.£ or $c.£j 

\\L^-^'-''-'[F,M,{^{L-\.)F2'){U]\\ < C2\\F[U (8.233) 

Proof. Since L, and ^„ are all commuting operators, any functions of these operators, defined by 

the spectral theorem, also commute. The spherical harmonic decomposition is into orthogonal subspaces 
preserved by these operators. Therefore, it is sufficient to prove the first two results for functions with a 
single spherical harmonic component and to prove the third result by considering the operator on the right 
as an operator on a single spherical harmonic. 
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For a fixed value of I, 

^'(r^ar)=*l(ZVl)*2(k|) (8.234) 

>x([l,oo),rVl)x([0,l],a;) (8.235) 

>x([r^oo),|a;|)x([0,l],a;) (8.236) 

>x{[l-\l],\x\) (8.237) 

>*i-^<|x|<i (8.238) 

Therefore, by the spectral theorem, for any function w e §, on each spherical harmonic 

ll*i-^<l€„l<i^ll < (8.239) 

Now the case when F2 is either $a,e or $c,£ is considered. Both ^a,e{x) and ^c,e{x) are smooth, 
strictly positive functions, so in either case ^2(2;) has a bounded inverse for \x\ G supp(\l/2) C [0,2], and 
(^2^2~^)(^n) is a bounded operator. Since all the operators involved commute and ^'i(Z''|a:|) is bounded, 
(\t'(L~^, •)F2~^)(^„) is a well defined, bounded operator. Again, for a fixed I, 

\x\>^l-'x{[ll-',oo),\x\) (8.240) 

>ll-'x{[l,^)J'\x\) (8.241) 

>i[-'^*i(/Vl) (8.242) 
>Cl-^'^i{l^\x\)(-^2Fi'){x) (8.243) 
>Cl-\'^{L-\»)F2^){x) (8.244) 

On each spherical harmonic, by the spectral theorem, 

ll^n^ll = IKI^nD^II > C\\L-\^{L-',.)F,-'){Uv\\ (8.245) 
Finally the commutator is calculated. Again this is proven on a single spherical harmonic shell. 

{^iL-',,)F,-'){U =*i(i'^J(*2F2"')(^n) (8.246) 

L'-^-''-'[F,{^^),{^{L-',.)F^-'){U]=L'-'^-''''[F,{^rn),^ (8.247) 

=i'-'"-"-'[i^i(x™),*i(L^CJ](*2F2-i)(^„) (8.248) 

+ i^-'"-"-'*i(L'x„)[Fi(ej, (*2F2"')(^ J] (8.249) 

Each of these terms is now estimated. Since only one spherical harmonic is being considered, the operator L 
can be replaced with the constant I. This simplifies the discussion of the commutator involving ^i{L^^^). 

Since ^'i is C°° and constant outside a compact interval, it follows that ^'^ is Schwartz class and that 
ll-^[*i]||i<oo. 

ll^'-'"-"-'[i^l(x,n),*l(/'U](*2i^2-')(^n)ll <Cl-' \\l'-'-[FM ,^ ,{1^ (8.250) 

<Cl-'\\Fi\U\ni^,{l'\.\)y]h (8.251) 

At this point, the derivative of the function in a scaled variable, {'i'i{l^»)y , is evaluated to be a scaled version 
of the derivative evaluated at the scaled variable, ^'^^'^ (/**•). To evaluate the Fourier transform of this, it is 
noted that ||^[/(A.)]||i = II J-[/(.)] 111. 

\\L'-^-''-'{F,{^m),^iiL'UmF2'){U\\ <Cr^||i^{||o.|K^.F[*i(/^|.|)]||i (8.252) 

<C||Fi'||^||.F[vl/l]||i (8.253) 

<C||i^i'||oo (8.254) 
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This completes the estimate on the first of the commutator terms. 

Since '^2 is Schwartz class, and F2 is smooth and has bounded inverse on supp(\E'2), (^'2^2^^) Schwartz 
class, and \\J^[i'i'2F2^y]\\i < 00. 

WL'-^^-'^.iL'UiFiM, {^2F2-'){U)\\ <C\\L^-"'--'[FM, {^2F2-'){U]\\ (8.255) 

<C||ii— "[Fi(x„0, (*2i^2"')(ejll (8.256) 

<C\\Fi\U\:F[i^2F2-'y]\\i (8.257) 

<C\\Fi\\^ (8.258) 

This completes the estimate on the second of the commutator terms. □ 

We now prove a l-^<\*\<i localised estimate. This is analogous to the $|»|<i localised estimate in 
lemma 15. 161 although, the proof is more complicated because the ^E" type localisations are more complicated 
to work with. Because ^„ is bounded below on the support of 5'l-«<|| |<i, we can dominate more powers 
of L with this localisation. Here, we take m = n. 



Lemma 8.19. For u ^ S is a solution to the wave equation, il + Hu = 0, e > 0, (5 > 0, and n G [0, i], 

(u, [H, r„,„]u) >\\L^-i-'-^L-^^\U<iXM? - 0{\\L-^x<.u\\\ L') (8.259) 



Proof. We take the result of lemmas [8. 151 and 18. 141 with m = n, 

d d 

(8.260) 

The term involving the derivatives is estimated first. The first step is weakening the estimate to eliminate 
the square root. 

\\L-^X^ixn)-^—^aAUXau\\ >\\U'-iX^ij,^) — <i>aAUXc.u\\ (8.261) 
op* op* 

Now is treated as if it were a localisation in and is commuted to the right of the chain of operators. 

||L^Xx(x„)— <J>a,.(|„)Xa^|| >||L"-t— Xx(x„)<I>a,e(OXa^^|| (8.262) 

- ||L"-t[Xx(x„),^]$,,,(^JXa«ll (8.263) 

The commutator from moving through Xj(x„) can be computed explicitly. It involves the function 
X'^{x) = —2xXi{x)'^, which is bounded. 

= -2L"x„Xx(x„)2 (8.265) 

=V'B. (8.266) 

The "localisation" L"--^ can be replaced by Li-'^(*(L"^ •)$-i)(|„) using lemma IHT^ 

||L^Xx(x„)^$,,,(^Jx.u|| >C\\L'-^-'i^{L-',.)^-^l){UXiM^aAUXc.u\\ (8.267) 

- WL^^'^iBxcuW (8.268) 

The (^'(L"'', •)$"),) (^„) is commuted back through the x„ localisation. The commutator involving 
(\E'(L~'', •)$~];)(^„) is shown to be lower order by lemma Lemma 18.131 can not be directly applied 
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to commutators involving {^{L •)^a,e)(^n)i because ^'(L ^,^n) ^ *i(-^'^^ri)^2(^„) is a product of two 
operators with different localisations. 

||L^lx(x„)^<I'a.e(l„)Xc.w|| >C\\L'-i-'Xi{^n)^iL-',UXc.u\\ (8.269) 

- \\L'-i-'[i^iL-', •)$-j)($J, li(x„)]<i>,,,(C JXaull (8.270) 

- ||i2"-4Bxaw|| (8.271) 
>C\\L'-i-'X^i^^)^SJiL~',UXo.u\\ (8.272) 

- C\\L^-i-^L^''-^+^Bxcu\\ - ||L2"^*Bxau|l (8.273) 

Since 2n — | is the power of L appearing in both error terms, and this exponent is less than ^^2 " ' ^^"^ error 
terms are 0{\\L~^ Xau\\^ , L^). To simplify the summation of the current estimate with the next step in 
this proof, an additional localisation by i—V[p~^) is introduced. 

||L^Xx(x„)A<i, (I (8.274) 
op* 

>C\\L'-^-'X^{^^){~VLp:')i^{L-',UXc.u\\ - 0{\\l'^Xo.u\\^L') (8.275) 

Now the gL"V[ terms are estimated. The goal is again to replace the localisation ^a.e[$,n) by ^[L^^ , ^„). 
This requires commuting {^[L^^ , •)^a,£)(^ii) through the p* localisation. A factor of is dropped to control 
commutator terms involving 5'(i^'',^,J at a later stage. 

\\L-^{-VipZ^)'-X^{^^)^aAUXcu\\ =||Li-t-^Xt(x„)(-y/p;i)^$^,,(|JXau|| (8.276) 
The bounded localisation (^(L~^, •)^~a\){^>n) introduced, and then commuted through the localisation in 



\\L-^{-Vip-'f^X^{^r,)<^aAUXo.u\\ (8.277) 

>C||(vI/(L-^.)$-i)(|„)Li-t-*XT(x„)(-V^/p-i)i$„,,(|JXa"ll (8.278) 

>Ci||Li-i-^XT(x„)(v|/(L-^.)$-i)(|J(-V^/p;i)^<i>,,,(^„)Xa^|| (8.279) 

-C72||Li-*-^[(f(L-^.)$-i)(^„),XT(x„)](-V^Zpr^)^$a,.(IJXa^l! (8.280) 

Now the operator (^'(L"'', •)$^];)(^„) is commuted through the localisation (— V^p^^)- This eliminates the 
<i>a,e(|„) localisation. 

\\L-^{^VipA)"-X^{^A^aAUXcu\\ (8.281) 

>C\\L^-i-'X^{^A{-VLpA)'^-^{L-\UXc.u\\ (8.282) 

-C||Ll-t-^XT(x„)[(*(L-^.)<i>-^)(^J,(-F/p;^)^]<i>a,.(l„)Xo^i|| (8.283) 

-||Ll-^-^[(*(i-^.)a>-l)(|„),XT(x„)](-Vlp:l)^ci>,,,(^Jx.^^ll (8.284) 

The commutator terms can now be estimated using lemma 15 . 1 81 and are found to be lower order error terms. 
Note that Xi^{x) has a i°° derivative. 

\\L^{~VipA)"^X^{^A'^aAtn)Xo.n\\ (8.285) 

>C\\L^-'^-'X^(^A(~VipA)--'i>{L-\Ux<.u\\ (8.286) 

- C||ii-i-^X|(x„)i"+^-iBi$a,,(ejxaw|| (8.287) 

- \\L^-i-'L^-+'-'BA-VipA)"^^aAUXc.u\\ (8.288) 
>C\\L^-^^-'X^{^A{~Vip-At^^{L-',UXc.u\\ (8.289) 

+ 0(||X — Xawf ,Li) (8.290) 
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Equations H8.275I) and 18.29UI) can be combined with the initial estimate on the expectation value of 
[H,Tn.n\ to produce an intermediate result. 

(u,[i7,r„,„]7i) >Ci\\L'-i-'X^i^n)i-V[p:')HiL-',UXc.uf (8.291) 
+ \\L'-^-'X^i^n)i-V[p:')H{L-',UXc.uf) (8.292) 
-Oi\\L^Xo.uf,L^) (8.293) 

Both these terms involve localisation in x„ acting on {—Vlp~^)i'if{L~^,^„)xaU. The sum of the locali- 
sations is bounded below by a constant. 



X.i^r + X,{xf = (^)^ + {x^J '-^r > C (8.294) 
Therefore, the two terms can be combined to provide a better estimate. 

(7.,[i?,r„,„]u) >c\\l'-^-\-vIp:')H{l-',Uxo.u\\ - oiWL'^xcuW^L') (8.295) 

To eliminate the factor of {—Vlp~^)^, a new function is introduced. This function, /, is a smooth, 
compactly supported function and equal to the inverse of {—Vlp^^)^ on suppxa. Since it is bounded, it can 
be freely introduced into the norm. The p^ localisation can then be commuted through 5'(L^'',^„). 

{u, [H,r„,^]u) >C\\L'-^-'f{-Vlp-')H{L-',UXau\\ - 0(||L^Xauf ,i') (8.296) 

>C\\L'-^-'^{L-',Ufi~Vlp:')h^u\\ (8.297) 

- \\L'-i-'[fi-V[p:')K^{L-',U]Xc.u\\ (8.298) 

-Oi\\L^X»u\\^,L^) (8.299) 

From the definition of /, the product of all the localisation reduces to Xa- The new commutator terms 
can be estimated by lemma lS.18l and are found to be lower order error terms. 

{u, [H, r„,„]7.) >C\\L'-i-'^{L-\ UXauW - \\L^-i-'L^+'-^Bxau\\ (8.300) 

-Oi\\L^Xc.u\\^L') (8.301) 

>C\\L'-^^~'^{L-\UXo.u\\ - OiWL'^XcuW^L^) (8.302) 

Lemma [8. 181 can now be used to replace the localisation in 5'(i^'',£„) by localisation in ^'i-i<|^ |<i. 

(u,[iJ,r„,„]u) >C||ii-4-^*i-.<|^j<iXau|| -0(||L^Xawf (8.303) 

□ 



8.5 Phase Space Induction 

The previous section shows that r„ i and r„_„ majorate L 2 and L 2 respectively. We would like to 
integrate the Heisenberg identity to conclude that the time integral of the expectation value of these powers 
of L are bounded. However, our definition of majoration allows the domination to occur only in a region of 
phase space and the lower order terms to be unlocalised. 

The lower order terms are 0{\\L~2— XauW^ , L^). To control these terms, we use a finite induction on n, 
to eventually control L^~^ without phase space localisation. 

Following this, we use the control of L^~^ in the conformal estimate to prove point wise in time, weighted 

decay, with bounds involving an additional factor. 

Theorem 8.20 (Phase Space Induction). If e > 0, then for u E E> a solution to the wave equation, 
il + Hu ~ 0, there is a constant, Kps, such that 

/oo 
\\L^-ixMtWdt < Kps{E[u,u] + \\u{l)f) (8.304) 
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Proof. This is proven for s = S + e. 

We induct on n to prove, simultaneously, the two statements 



Ili^^^Xauf =0{L'), (8.305) 

||ii-V[i'"",oo), ^hcuf ^0{L'). (8.306) 

The induction will run from n = 0, in steps of size 5, as long as 

n < i - 2(5 - e. (8.307) 

Each step in the induction will be proven by Morawetz type arguments using r„^„ and r„ i . 

The base case of (|8.305f) follows from the angular modulation result, theorem 17.61 which says, in the 
notation of this section = 0{L^). The base case of H8.306|l follows from the angular modulation 

theorem, theorem 17.61 For n = 0, since 

l^-'x{[l, oo), x) <l^-'x{[l^'', oo), x) < X, (8.308) 
by the spectral theorem and corollarv l7.7l 

\\L'-'x{[L,^), 4-)Xcur <||/-Xau|p = 0{L^). (8.309) 

The inductive step is now considered. 

Morawetz type estimates with r„_„ and r„ i will be used. Before these estimates are proven, ||r„^,„M||^ 
is shown to be bounded by the energy and a local decay term, under the condition < n < m < i. 

To begin estimating the norm of r„ „jM, r„ „ is expanded as a product of operators, and the factor of 
7iTO is expanded as a sum of two terms. These terms can be further simplified by eliminating bounded 
functions. 

||r„.,„«|| =||Xa<fa,e(C„)7L'"i"~'$a,e(eJXaU|| (8.310) 
<||Xa<i>a,e(^jT^.9L".i"-^<i>a,.(l„)XoU|| (8.311) 

op* 



1 L" 

xAAUt; , ,.L^ L"-'^aAUXo.u\\ (8.312) 
d 



<ll$a,.(l„)^5L".i"-^<fa,.(eJXaU|| + ||£™+"-^Xa^|l (8.313) 

Interpolation can be used to simplify \\L''^^^^'^Xau\\. 

||r„,™u|| <|j$^^^(|J_g^,„L»~^$,_,(|J;^„^|| + C{E[u] + II (1 + piy' u\\') (8.314) 

To control ||$a,£(^„)g^5L'"i""'^*i'a,e(^„)Xaw||; it IS rewritten in terms of ^„ and L, and its sub factors are 
rearranged. 

\\<^aAU-g^^L"-'9L^^aAUXc.u\\ =\\L'~'<faAU^n9L'^'^aAUXaU\\ (8.315) 

<\\L'-'9L^^^^n^aAU^XaU\\ (8.316) 

+ \\L'-'[^„^aAU,9LA^aAUxJ (8.317) 

The first of these two norms can be estimated by lemma IHTtI The second can be estimated using lemma IH.13I 
which states that the commutator is L'^+^^^B. 

\\<^aAUj^L^~"9L^'^aAUXc.u\\ <\\L-^^Xo.u\\ + \\j—Xo.u\\ (8.318) 
ap* op* 

+ ||L"+™-^i?<&,,,(eJxau|| (8.319) 
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Since "'^"'""1^ < 1 and n + m — e < 1, each of the terms involving powers of L in this expression can be 
estimated by interpolation. 



l^aAU^^L^-'gL'^^aAUXaUW <CiE[u] + \\ {l + pi) ' uf) 



This completes the estimate on ||r„_„iu|j. 

We can now integrate the Heisenberg like identity to find 



{u, [H,Tn,m]u)dt 



<C{E[u] + II (f + piy' u{T)r + II (1 + pI) u{l)f) 



(8.320) 



(8.321) 

(8.322) 
(8.323) 



The second term on the right can be dropped since, by the local decay result, the norm || (l + pj) u{T)\\^ 
on a sequence of times. 



{u,[H,r.n,,n]u)dt<C{E[u] + \\ {I + pD 'm(I)IP) 



(8.324) 



The derivative localisation results, lemmas 18.191 and 18.161 can now be applied. From the inductive hy- 
pothesis (|8.305|l and condition H8.307|) . the lower order terms are integrable: 



\\L'^Xc.ur=0{L'), 
0{\\L^Xc.u\\\L^) =0{L^). 
From inductive hypothesis (|8.3U6|I and the second derivative localising result, lemma IS. 191 

||Li-^x([i'-"-',oo),^)Xau|P=||i^-'x([i'-"-',i'-"),7^)Xow||' 

+ ||Li-^x([i'-",oo),^)Xau||2 

= ||il-*'*L-.<|^J<lXaM|P 



9p* 

< {u, [H, r„,,» + 0(||L^Xawf , L') 
+ 0{L^). 



Integration in time extends inductive hypothesis IjS. 306(1 to n + 6 



d 



\\L'-'x{[L^^'''^\^).i-)Xo.ufdt<C{E[u] + \\u{l)f). 
op* 



(8.325) 
(8.326) 

(8.327) 

(8.328) 
(8.329) 
(8.330) 

(8.331) 
(8.332) 

(8.333) 



From condition (|8.307|l . it follows that '^^'^ < \ — 6. From the first derivative localising result, lemma 



"Xa^ilP =||L^X([0,Ll-"], A)^^y||2 

dp* 

+ ||L^x([i'-",oo),^)Xa«ir 
ap* 



<||L^$I« j<iXau|r + ||Li-'x([i'-", oo), -^)Xc.nf 
<{u, [H,r^i]u) + 0{\\L'^Xc.nf,L') + 0{L'). 



(8.334) 

(8.335) 

(8.336) 
(8.337) 
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Integration in time extends inductive hypothesis H8.3U5|) up to n + S. 

||L^T— XaulPdt <C{E[u] + ||u(l)||2) (8.338) 



Since the induction continues as long as condition l|8.307|) holds, after the last application of equation 



\\L — '^Xaurdt <C(E[u] + ll'u(l)lr) (8.339) 

This proves the desired result with e = 5 — e. Since 5 and e can be taken to be arbitrarily small, so can e. □ 

Except for the loss of L'^, this provides the control on the angular energy near the photon sphere required 
by corollarv l7.8l to control the weighted norm. Since the wave equation is linear and L commutes with 
H , L'^u is also a solution, and the energy of this function can be used to recover the additional factor of L'^u. 

Theorem 8.21. If u is a solution to the wave equation H2.21|l , and u — ru, then 

II J^^u|Il6(ot) <t'^'C[\\uo\\l + E[uo, ui] + Ec[uo, ui] + \\L'uf + E[L'uo, L'ui])i (8.340) 

Wipl + l)^u\\LHm) <t'^Ci\\uo\\l + E[uo, ui] + Ec[uo, ui] + \\L'u\\^ + E[L'uo, L'ui])^ (8.341) 

Proof. If M is a solution to equation (|2.21|) , then u is a solution to equation (|2.24|l . The weighted norms 
are related. 



- J F^\ufd^n (8.342) 
: J F^r~^\u\d^n. (8.343) 



Since the wave equation equation H2.24|l is linear and H commutes with L, if w is a solution, then so is 
L'^u. The phase space induction theorem, theorem 18. 201 can be applied to L'^u. 

\\L^-^Xc{L'u)r =0{L^), (8.344) 
J \\Lxau\\^dt <C(||i'=u||2 + E[L'u, L'ii]). (8.345) 

This is the term which is needed to control the conformal charge in the reduction to angular energy on 
the photon sphere, corollarv l7.8l 

Ec[u{t),u{t)] <Ec[uo,ui] + CtE[uo,ui]^E[uQ,ui]i + \\uo\\) (8.346) 

+ CtJ j Xa\^s2u\'^d^ijLdT (8.347) 

<Ec[uQ, Ml] + CtE[uo, u{\^{E[ua, ui]^ + lluoll) (8.348) 
+ Ct{\\L''uof + ElL^uo, L^ui]) (8.349) 
Since the conformal charge controls the weighted norm, we have 

WF^^^Wmm) <( / i^V-4|M|dV)^ (8.350) 
<Ct'^{\\uof + E[uo,tii]+Ec[uo,ui] + \\U-uaf +E[L''uo,U-ui]). (8.351) 
Similarly, since lemma controls the weighted norm by the conformal charge, 

-^^r^'^)L^m) ^^J^") (8.352) 
<Ct-'^Ec[u{t),u{t)] (8.353) 
<Ct-^Ec[uo, ui] + Ct~i(i;[uo, Ml] + IImoII^ + Wuof + S[L^uo, i'ui]) (8.354) 

□ 
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